AO-A064  944  GENERAL  ELECTRIC  CO  CINCINNATI  OHIO  AIRCRAFT  ENGINE  GROUP  F/G  21/5 


ANALYTICAL  DERIVATIVES. (U) 


ABA06494 


FINAL  TECHNICAL  REPORT 


General  Electric  Company 
Cincinnati,  Ohio  45215 


December  1978 


AIR  PORCE  AERO  PROPULSION  LAtORATORY 
AIR  PORCI  SYSTEMS  COMMANO 
WRIOHT- PATTERSON  APB.  OHIO  4S433 


w 


NOTICE 


When  Government  drawings,  specifications , or  other  data  are  used  for  any  pur- 
pose other  than  in  connection  with  a definitely  related  Government  procurement 
operation,  the  United  States  Government  thereby  incurs  no  responsibility  nor  any 
obligation  whatsoever;  and  the  fact  that  the  government  may  have  formulated, 
furnished,  or  in  any  way  supplied  the  said  drawings , specifications , or  other 
data,  is  not  to  be  regarded  by  implication  or  otherwise  as  in  any  manner  licen- 
sing the  holder  or  any  other  person  or  corporation,  or  conveying  any  rights  or 
permission  to  manufacture,  use,  or  sell  any  patented  invention  that  may  in  any 
way  be  related  thereto. 


This  report  has  been  reviewed  by  the  Information  Office  (OI)  and  is  releasable 
to  the  National  Technical  Information  Service  (NTIS) . At  NTIS,  it  will  be  avail- 
able to  the  general  public,  including  foreign  nations. 


This  technical  report  has  been  reviewed  and  is  approved  for  publication. 


' i s 

JAMES  R.  RUBLE 

Technical  Manager  JACK  RICHENS,  Chief 

Turbine  Engine  Performance  Performance  Branch 


FOR  THE  COMMANDER 


ERNEST  C.  SIMPSON,  Director 
Turbine  Engine  Division 


"If  your  address  has  changed,  if  you  wish  to  be  removed  from  our  mailing  list, 
or  if  the  addressee  is  no  longer  employed  by  your  organization  please  notify 
AFAPL/TBA  ,W-pafb,  OH  45 433  to  help  us  maintain  a current  mailing  list". 


Copies  of  this  report  should  not  be  returned  unless  return  is  required  by  se- 
curity considerations , contractual  obligations , or  notice  on  a specific  document. 

AIR  FORCE/56780/9  February  1979  — 55 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  of  THIS  PAGE  (Wb*n  DM*  Enttttd) 


lEPORT  DOCUMENTATION  PAGE 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FOR.-4 


MergrtT>«&M«^ -7^ 

/AFAPL-|TR-78-j/l  X 


[1.  OOVT  ACCESSION  NO 


S RECIPIENT'S  CATALOG  NUMBER 


4.  TITLE  (mid  Submit) 


ANALYTICAL  DERIVATIVES. 


i>  fTMPB  BF  RBPBRP  4 PBRIB8  CBMBRSB 

Final  /echnical  - 


l /epmn 


r:  * 

D.F./Berg  \ 
W.C. /Colley  \ 
erse  \ 


R78AEG517/ 

"WWWTCT'BR’grant  NUMiE*r<j 


F33615-78-C-22J03  «Vt &Of 


G.L 


.JO 


onverse 


W 


16.  PROGRAM  EL  EMC  NT.  PROJECT.  TASK 
AREA  A WORK  UNIT  NUMBER* 

^Jj^.ll-35  (Pj 


« ORGANISATION  NAME  ANO  ADDRESS 


General  Electric  Company 
Aircraft  Engine  Group 
Cincinnati,  Ohio  45215 


I I.  CONTROLLING  OFFICE  NAME  ANO  ADDRESS 

Air  Force  Aero  Propulsion  Laboratory  (TBA)  ^ / / 
Wright-Patterson  AFB,  Ohio  45433 


G2 


De 


78 


" nnnn III ~ "Til— r n *I1T  * i~ITf  tlflirr  - If  t a—*'*--  ■ 


IS.  SECURITY  CLASS.  (•!  tbit  r**onj 

Unclassified 


IS*,  oeclassification/downgrading 
^M^Le  N/A 


16.  DISTRIBUTION  STATEMENT  ( o I Ih It  Rtporl) 

Approved  for  public  release:  distribution  unlimited 


17.  DISTRIBUTION  STATEMENT  fofjfi*  abetted  entered  In  Block  20,  II  dlllerent  from  Report) 


D D C 


lIHfrOlHHU-pm 

FEB  27  1979 


is.  supplementary  notes 


ilitolLU  FB 

B 


19.  KEY  WORDS  (Continue  on  rtmit  elde  It  neceeeery  end  Identity  by  block  number) 

Analytical  Derivatives 
Engine  Cycle  Balance 
Iteration  Techniques 


|S0.  ABSTRACT  (Continue  on  revet  ee  elde  If  neceeemry  end  Identity  by  block  number) 

--The  variable  cycle  engines  now  being  studied  for  preliminary  engine-aircraft| 
design  have  greatly  increased  the  complexity  of  engine  modeling,  control 
element  design,  and  cycle  calculation  relative  to  the  simpler  turbofan  and 
turbojet  engines  which  preceeded  them.  The  production  of  engine  data  is  the! 
major  computations  cost  element  in  mission  studies  and  optimization  investi-f 
gations;  this  problem  can  become  even  more  severe  in  the  future  unless  1 — >, 


EDITION  OF  • NOV  6>  I*  OBSOLETE 


SECURITY  CLASSIFICATION  OF  THIS  PAOE  (Mhtti  Dm*  B*l*t*d) 


<L<~ 


DO  ,:srw  1473 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  This  F AGKPYhan  Dmf  Ini. fMO 


V) 

improvements  in  cycle  calculations  and  component  modeling  techniques  are 
developed  and  introduced  now. 

The  current  method  of  obtaining  balanced  cycle  data  points  is  to  use  the 
simultaneous  Newton-Raphson  iteration  method.  The  partial  derivatives  are 
obtained  by  numerical  differentiation.  This  report  documents  the  results 
obtained  by  calculating  the  partial  derivatives  from  analytical  expressions 
obtained  by  differentiating  a cycle  deck.  Examples  illustrating  both  the 
method  of  obtaining  the  analytical  derivatives  as  well  as  setting  up  the 
control  logic  are  given.  A cost  comparison  was  carried  out  between  two 
engine  simulations  using  numerical  derivatives  and  the  same  pair  of  engine 
simulations  using  analytical  derivatives.  A data  matrix  consisting  of  the 
same  411  operating  points  was  used  for  each  of  the  engine  simulation 
comparisons.  ^ 

A cost  saving  of  about  44  percent  was  obtained  for  the  deck  run  internally 
and  52  percent  for  the  deck  run  externally  (the  difference  is  due  to  the 
use  of  a larger  set  of  output  parameters  for  the  internal  deck) . It  is 
concluded  that  the  greatest  return  per  dollar  of  cost  results  when  only 
that  portion  of  the  model  for  which  the  coding  remains  relatively  constant 
j is  differentiated.  Applying  this  philosophy  an  annual  saving  of  about 
$75,000  is  estimated  for  the  running  of  internal  (AEG-Evendale)  cycle  decks. 


UNCLASSIFIED 


ICCUR'TV  CLASSIFICATION  OF  THIS  FAOKr»Mn  Of  lMm« 


PREFACE 


This  report  describes  a design  study  effort  conducted  by  the  General  Electric 
Company  and  sponsored  by  the  Turbine  Engine  Division  of  the  Air  Force  Aero 
Propulsion  Laboratory,  Air  Force  Systems  Command,  Wright-Patterson  AFB,  Ohio 
under  Project  3066-11-35,  Contract  F33615-78-C-2203.  Mr.  James  R.  Ruble,  AFAPL/ 
TEA,  va6  the  Air  Force  Project  Engineer.  

The  work  reported  herein  was  performed  during  the  four-month  period 
beginning  June  1978  and  ending  September  1978.  The  GE  Engineering  Manager  was 
Donald  F.  Berg  who  was  assisted  principally  by  William  C.  Colley  and  George  L. 
Converse.  The  Program  Manager  was  Donald  E.  Uehling. 
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SUMMARY 


The  next  generation  of  military  aircraft  weapon  systems  will  be  required 
to  achieve  substantial  improvements  in  design  mission  performance  and  effec- 
tiveness, multimission  versatility,  life-cycle  costs,  and  survivability. 

These  requirements  will  place  increased  demands  on  the  propulsion  system  for 
advanced  cycle  concepts,  advanced  material  and  design  technology,  variable 
geometry  capabilities,  and  a more  effective  engine-airframe  installation  with 
particular  emphasis  on  inlet-engine  airflow  matching  and  engine-airframe 
thrust  matching  across  the  complete  operating  regime.  The  resulting  propul- 
sion system  concept  and  evaluation  and  cycle  selection  process  will  require 
a more  complex  engine-airframe  interaction  which  will  involve  a substantially 
greater  number  of  engine  and  airframe  design  parameters,  more  extensive  con- 
trol and  scheduling  requirements,  and  a broader  spectrum  of  mission  and  oper- 
ational requirements. 

The  studies  necessary  to  satisfy  these  increased  design  requirements  will 
necessitate  calculating  many  more  balanced  engine  cycle  performance  points. 
This  study  program,  Analytical  Derivatives,  was  aimed  at  reducing  the  computer 
cost  for  calculating  engine  cycle  performance  by  utilizing  analytical  deriva- 
tives instead  of  finite  difference  derivatives  in  the  engine  cycle  balance 
iteration  procedure.  The  technique  developed  is  to  differentiate  the  engine 
cycle  and  to  use  the  numerical  values  of  these  analytical  derivatives  in  the 
cycle  balance  iteration  procedure.  The  calculation  of  the  analytical  deriva- 
tives will  require  less  computer  processor  time  than  the  calculation  of  fi- 
nite difference  derivatives  which  require  evaluating  the  complete  cycle  for 
each  independent  variable.  This  technique  was  applied  to  a Variable  Cycle 
Engine  simulation  and  the  resulting  computer  program  was  used  to  calculate 
411  flight  operating  points.  The  results  indicated  that  the  computer  costs 
were  reduced  about  52Z  when  using  the  WPAFB  CDC6600/ CYBER  74  computer. 

During  this  study,  an  alternate  approach  of  applying  the  analytical 
derivative  concept  was  devised.  The  differentiation  of  the  complete  engine 
cycle  requires  the  constant  updating  of  these  analytical  derivatives  whenever 
a new  component  is  added  or  equations  are  changed  or  added  to  an  existing 
component.  This  alternate  approach  differentiates  only  the  low  level  sub- 
routines, such  as  those  subroutines  which  calculate  thermodynamic  and  aero- 
dynamic properties.  These  subroutines  represent  approximately  10Z  of  engine 
cycle  program  logic  but  account  for  about  75Z  of  the  computer  processor 
time.  By  applying  the  analytical  derivative  technique  the  calculation  time 
for  these  subroutines  will  be  reduced  and  since  these  subroutines  remain 
unchanged  changes  in  the  engine  cycle  will  not  require  additions  to  this 
set  of  analytical  derivatives.  Utilizing  the  same  Variable  Engine  Cycle 
program  to  calculate  the  411  flight  operating  points,  the  results  showed 
about  a 47Z  reduction  in  computer  processor  time  on  a Honeywell  6000  com- 
puter. Using  the  ratio  of  the  Honeywell  6000  computer  time  savings  for  the 
alternate  method  to  the  analytical  derivative  technique,  it  is  estimated 
that  the  alternate  method  would  result  in  a 40Z  costs  savings  using  the  WPAFB 
CDC6600/CYBER74  computer. 


It  can  be  concluded  that  in  any  engine  simulation  model  using  numerical 
derivatives  to  obtain  cycle  points,  a considerable  savings  in  both  computer 
time  and  cost  can  be  obtained  by  the  inclusion  of  analytic  derivatives.  The 
greatest  return  per  dollar  of  cost  results  when  only  that  portion  of  the 
model  for  which  coding  remains  relatively  constant  is  differentiated,  ^or 
the  Air  Force  engine  simulation  model  (SMOTE),  the  entire  deck  should  be  dif- 
ferentiated since  the  coding  changes  do  not  appear  to  occur  frequently. 
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1.0  INTRODUCTION 


The  Aircraft  Engine  Group  (AEG)  of  the  General  Electric  Company  has  com- 
pleted a contract  entitled  "Analytical  Derivatives"  which  was  sponsored  by  the 
Turbine  Engine  Divison  of  the  Air  Force  Aero  Propulsion  Laboratory.  The  Air 
Force  Project  Engineer  was  James  R.  Ruble  of  the  AFAPL  Performance  Branch. 

The  overall  objective  of  the  contract  effort  was  to  develop  and  assess  the 
suitability  and  the  cost  reduction  potential  of  the  analytical  derivative 
procedure  as  applied  to  the  computer  modeling  of  turbine  engine  performance. 

In  particular,  this  program  developed  the  relevant  mathematical  expressions 
for  the  analytic  derivatives  and  applied  the  analytic  derivative  technique  to 
the  calculation  of  balanced  cycle  performance  for  a selected  Variable  Cycle 
Engine  (VCE)  simulation. 
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2.0  SELECTION  OF  ENGINE  MODEL 


The  analytical  derivative  study  utilized  tvo  point  design  Variable  Cycle 
Engine  (VCE)  simulations  which  are  AEG  inhouse  computer  prograsis.  A schemat- 
ic of  the  VCE  used  during  the  study  is  shown  in  Figure  1.  Both  the  fan  blocks 
and  the  compressor  representation  within  the  VCE  simulation  exist  in  parametric 
form.  The  remainder  of  the  cycle  equations  are  consistent  with  those  utilized 
in  parametric  cycle  decks,  thereby  assuring  that  all  analytical  derivatives 
generated  in  the  program  effort  would  be  useful  in  full  parametric  performance 
decks.  By  utilizing  the  parametric  feature,  the  engine  can  be  changed  easily 
by  selecting  different  compressor  and  fan  designs  from  the  family  of  designs 
available  in  the  model. 

The  demonstration  of  the  analytical  derivative  procedure  was  carried  out 
using  twc  different  engines.  The  initial  evaluation  used  an  engine  with  a 
fan,  core-driven  booster,  and  a compressor,  each  selected  from  the  parametric 
system  by  design  pressure  ratios  of  2.97,  1.35,  and  7.00,  respectively.  This 
was  designated  as  Engine  1.  For  the  second  engine,  the  design  pressure  ratios 
of  these  components  were  3.28,  1.52,  and  5.60.  The  second  engine  was  desig- 
nated as  Engine  2. 


4 


42  49  15B 


3.0  DESCRIPTION  OF  THE  ANALYTICAL  DERIVATIVE  METHOD 


In  order  to  balance  an  engine  cycle,  it  is  necessary  to  solve  a set  of 
n simultaneous  non-linear  algebraic  equations  in  n independent  variables, 
represented  symbolically  by: 

^i(Xt , X2,  ...  Xn)  ■ 0 i - 1,  2,  ...  n 

Since  the  equations  are  nonlinear,  direct  solution  is  not  possible.  The 
solution  is  found  by  trial  and  error,  using  the  Newton-Raphson  method  to 
guide  convergence.  This  method  is  summarized  briefly. 

In  general,  for  some  set  of  approximate  values  of  the  independent  vari- 
ables, the  equations  will  not  be  satisfied;  instead,  the  values  of  the  func- 
tions Y{  will  differ  from  zero  by  some  error 

YiUi,  X2,  ...  Xn)  - Ei  i - 1,  2,  ...  n 

The  object  is  to  determine  a set  of  corrections  4X£  which  will  force  the 
erors  to  zero: 


Y. 

l 


♦ 


3Yi 

3X^ 


«X, 


«x„ 


♦ 


«x 

n 


0 i 


1.  2, 


• • • n 


The  Y^'s  are  evaluated  using  the  aproximate  values  of  the  X'a.  In  princi- 
pal, the  partial  derivatives  can  also  be  evaluated,  so  that  the  above  equa- 
tions form  a set  of  simultaneous  linear  equations  which  can  be  solved  di- 
rectly for  the  corrections  SX£.  The  effectiveness  of  the  corrections  depends 
upon  the  accuracy  of  the  approximations  to  the  root  values  of  the  X‘s  and  the 
linearity  of  the  original  simultaneous  equations. 

In  engine  cycles,  the  functions  Yi  are  so  complicated  that  the  partial 
derivatives  have  usually  been  evaluated  numerically: 


3Y*  _ YjUi,  X2,  ...  Xj  ♦ AXj,  ...  Xn)  - Y j(Xj , X2,  ...  X j , ...  Xn) 
TxJ  5 (Xj  ♦ AXj)  - Xj 

This  procedure  requires  n ♦ 1 complete  passes  through  the  engine  model  to 
compute  each  set  of  Y’s  and  partial  derivatives. 
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Direct  evaluation  of  the  partial  derivatives  requires  much  less  compu- 
tion  time,  as  will  be  shown,  although  a completely  different  logic  path 
through  the  engine  cycle  is  required. 


3.1  ANALYTICAL  DERIVATIVE  APPROACH  - COMPLETE  DERIVATIVE  SET 

A description  of  the  basic  analytical  derivative  approach  follows. 

Each  component  module  has  a list  of  needed  inputs  (Xj,  X2 Xn)  and  calculates 

a set  of  dependent  variables  (Y^,  Y2,  Y3 Ym)  utilized  in  other  modules. 

A set  of  total  differential  equations  can  be  mathematically  derived: 


dX 


dX 


dY 

„ ®Y1 

dX, 

f 3Y1 

dX„ 

-♦ill 

1 

axx 

1 

+ ax2 

2 

axn 

dY 

- 3Y2 

dX 

9Y2 

dX„ 

— ♦!£ 

2 

3XX 

1 

ax2 

2 

3Xn 

dY 

. 3Yn> 

dX 

dX„ 

m 

3XX 

1 

+ ax2 

2 

aXn 

The 

partial  derivatives 

aYx/aXi,  ; 

groupings  of  parameters  calculated  on  the  base  point.  The  total  differential 
equations  can  be  derived  and  programmed  in  each  subroutine.  When  the  sub- 
routines are  executed  in  order,  the  values  of  dY^  from  an  upstream  sub- 
routine frequently  become  the  input  values  of  dX^  needed  in  a downstream  cal- 
culation. In  this  way,  the  derivative  expressions  can  be  modularized  and  uti- 
lized in  a variety  of  cycles  consisting  of  different  ordered  sets  of  sub- 
program modules. 

To  iterate  a cycle  to  balance,  partial  derivatives  with  respect  to  a set 

of  independent  variables  , Z2,  Zn  are  required.  This  set  of  variables 

is  a sub-set  of  total  set  of  Xj's  utilized  in  each  subroutine  as  independent 
variables.  Similarly,  a subset  of  the  Y^'s  represents  a set  of  variables  to 
be  driven  to  zero  by  the  iterative  technique.  To  evaluate  the  partial  deriv- 
atives of  the  Y's  needed  by  the  iterative  technique  with  respect  to,  say,  Zj , 
the  value  of  dXj  corresponding  to  X{  is  set  equal  to  unity  and  all  other  dX 
values  corresponding  to  Z2  through  Zn  are  set  to  zero.  By  calculation  through 
the  derivative  equations  the  appropriate  derivative  values  can  be  evaluated. 

The  derivative  equations  must  be  evaluated  in  this  manner  once  for  each  inde- 
pendent iteration  variable.  Although  this  procedure  leads  to  multiple  evalua- 
tions of  the  derivative  equations,  it  does  make  the  derivative  equations  inde- 
pendent of  subroutine  order  and  only  dependent  on  a fixed  set  of  X^'s  and 
Yj's  in  each  module  (unique  interface  parameter  sets  between  routines).  More 
importantly,  perhaps,  it  permits  the  derivative  equations  to  be  programmed  in 
a much  more  compact  manner.  Large  portions  of  derivative  evaluations  are 
common  to  different  independent  variables  and  need  not  be  recalculated  for  each 
variable . 
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As  an  example,  consider  a core  compressor.  Two  of  the  input  variables 
(X's)  could  be  core  speed  (PCN2S)  and  inlet  air  temperature  (T25).  Corre- 
sponding output  variables  (Y's)  would  be  stall  margin  (SM25)  and  outlet  air 
temperature  (T3).  The  core  speed  is  typical  of  the  independent  iteration 
variables  (Z's)  while  T25  is  calculated  in  an  upstream  component.  Similarly, 
stall  margin  is  typically  a dependent  iteration  variable  (driven  to  a demand 
value)  while  T3  is  used  as  an  input  to  the  combustor.  The  equations  would 
take  the  form: 


dT3 


dSM25 


3T3 
3T25 

3SM25 

3T25 


dT25  ♦ 


3T3 


dT25  + 


3PCN25 
3SM25 


dPCN25  ♦ ... 


3PCN25 


dPCN25  ♦ ... 


With  dPCN25  set  equal  to  unity  and  the  differentials  of  all  other  independent 
iteration  variables  dZ's  set  equal  to  zero  the  set  of  equations  reduces  to: 

dT3  = 3T3/3  PCN25 

dSM25  = 3SM25/3  PCN25 

If,  on  the  other  hand,  we  set  dPCN25  equal  to  zero  and  set  the  differential 
of  some  other  independent  variable  (say  fan  speed  dPCN2)  to  unity,  then  the 
equations  reduce  to: 


dT3 


dSM25 


3T3 


3T25 

3SM25 

3T25 


dT25 


dT25 


Although  the  magnitude  of  dT3  is  different  for  the  two  cases,  the  same  com- 
bustor and  other  downstream  equations  (utilizing  dT3)  can  be  used. 

Note  that,  in  the  second  example,  the  temperature  T2S  is  a function  of 
fan  speed,  so  that: 

dT25  • -nib dpCH2  * ••• 

through  the  chaining  of  logic  in  preceeding  component  modules. 

1,  and  all  other  independent  variable  differentials  are  zero. 


When  dPCN2 


dT3 


3T3  3T25 


3T3 


3T25  3PCN2  3PCN2 
although  the  terms  are  not  actually  collected  in  this  form. 


8 


3.2  SUBROUTINE  DERIVATIVE  STRUCTURE 


in  general,  Che  toCal  differential  equations  described  previously  were 
programmed  into  the  high-level  (component  module)  subroutines,  even  though 
evaluation  of  some  variables  occurs  in  low-level  utility  routines.  The 
utility  rouitines  were  modified  as  required  to  provide  values  of  partial 
derivatives.  Exceptions  were  the  compressor  and  turbine  map  subroutines, 
where  the  differentiation  was  performed  in  the  subroutines. 

Differentiation  of  most  of  the  formulae  was  straightforward,  since, 
usually,  a dependent  variable  was  expressed  explicitly  as  a function  of 
independent  variables: 

Y - f(Xlt  X2,  ...  Xn) 


dY  ■ Hr  d*i  * ¥r2  ‘h  * 


dX 

n 


Even  in  cases  where  the  equations  could  not  be  solved  explicitly  for  the  de- 
pendent variables,  and  iteration  was  required  for  evaluation  of  the  variables 
themselves,  the  derivatives  of  the  functions  could  always  be  obtained,  so 
that  explicit  linear  relations  for  the  differentials  were  obtained,  which  did 
not  require  iteration. 

For  a system  of  equations: 

f i(Xi,  x2,  •••  Xn,  Yj,  Y2,  ...  Ym)  “0  i " 1,  2,  ...  m 


the  total  differential  expressions  reduce  to: 


3f;  3f i 

IyT  dYl  + dY2  + 


♦ 


3fi 

*Ym 


dY 

m 


dX, 


+ 


3fi 

3X2 


dX„ 


+ 


♦ 


i * 1,  2,  . . . m 


which  are  linear  in  the  dy's  and  can  be  solved  directly. 


The  algebraic  terms  representing  the  partial  derivatives  in  the  above 
formulae  are  evaluated  during  the  first  derivative  pass,  using  values  of 
independent  and  dependent  variables  obtained  during  the  base  pass.  Since, 
in  subsequent  derivative  passes,  only  the  dX's  change,  the  values  calculated 
in  the  first  pass  are  retained  in  storage  arrays  for  reuse.  Additional  com- 
putation time  was  saved  by  performing  derivative  calculations  in  only  part  of 
the  engine  model  on  each  derivative  pass.  Derivative  calculations  downstream 
of  the  point  where  the  last  dependent  iteration  variable  is  evaluated  are 
unnecessary,  as  are  calculations  upstream  of  the  point  where  each  independent 
iteration  variable  first  appears,  except  for  evaluation  of  partial  derivative 
terms  on  the  first  derivative  pass.  First-guess  and  engine  sizing  formulae 
were  not  differentiated. 


In  the  General  Electric  cycle  decks,  three  different  thermodynamic 
property  subroutines  are  used,  in  each  of  which  the  variables  by  InP  (log  of 
pressure),  FAR  (fuel  air  ratio)  and  WAR  (water  air  ratio)  are  independent. 

Of  the  variables  T,  H and  S,  one  is  independent  and  the  other  two  dependent, 
according  to  which  subroutine  is  called.  Other  dependent  variables  are  R (gas 
constant),  CS  (ratio  of  specific  heats  at  contant  entropy),  ALPHA  and  BETA 
(defined  in  the  following  equations).  In  differentiating  these  relations, 

WAR,  CS,  ALPHA,  and  BETA  were  considered  to  be  constants.  The  differentia- 
tion formulae  were  derived  as  follows: 


H =*  H (T,  InP,  FAR) 


dH 


ill 

3T 


dT  ♦ 
P,  FAR 


3H 

(HpT 


| d(lnP) 
|T,  FAR 


+ 


3H 

3 (FAR) 


d(FAR) 

T,P 


by  standard  thermodynamic  relations,  (Reference  1): 


3H\ 

3VP,  FAR 


CP  * CS  (1  - ALPHA) 


ALPHA  - 


FAR 


3H _ p/3H\ 

3 (InP)  T|  FAr  \aP/T,  FAR 


(ALPHA) (R)T 


The  derivative  functions  CS  and  ALPHA  were  already  returned  from  the  thermo 
routines.  The  coding  was  modified  to  compute  and  return  the  remaining  deriv- 
ative 3H/3(FAR)tjp- 


Thus , 

dH  * CS  (1  - ALPHA)  dT  + ( ALPHA) RT  d(lnP)  ♦ 


3H 

3 (FAR) 


T,P 


d (FAR) 


This  formula  is  readily  inverted  to  exchange  dependent  and  independent 
variables : 

dH  - ( ALPHA )RT  d(lnP)  - |jFAR)  d(FAR) 

dT  - , LlfP 

CS  (1  - ALPHA) 


Similarly: 

S - S (T,  InP,  FAR) 


dS  * || 

3T 


dT  + I7TSpV 
FAR.P  »ClnP) 


d(lnP)  ♦ 

FAR.T  3 (FAR)  Tp 


d(FAR) 


dS  - CS  (1  - ALPHA)  - R (1  - ALPHA)  d InP  ♦ g'^FAR)  d(FAR) 

T,P 
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and 


d (FAR) 
T,  P 


dR  - H 


dT 


3R 


3 (lnP) 


d(lnp)  + TTfarT 

T,  FAR 


|P, FAR 

using  the  definition  of  ALPHA  above  and 


BETA 


P m _ I \ 

R vp/t,  far  R \3(lnPy  x, 


FAR 


then 


dR  = - (ALPHA)R  ^ - (BETA)R  d(lnP)  + 


T,P 


d (FAR) 


For  lean  gas  mixtures  not  subject  to  chemnical  dissociation,  ALPHA  ** 
BETA  * 0,  and  the  thermodynamic  properties  at  a given  temperature  become 
linear  functions  of  fuel-air  ratio  when  expressed  per  pound  of  dry  air;  for 
example : 

H - Ha  + (FAR)  Hf  + (WAR)  H*, 

3H 


3 (FAR) 


- H 

T f 


where  Hf  is  the  enthalpy  of  the  CO2  and  H2O  formed  by  burning  a pound  of 
fuel,  less  the  oxygen  consumed  and  is  a function  of  temperature  only,  and  Hw 
is  the  enthalpy  of  the  water. 

Other  low-level  utility  subroutines  for  which  derivative  formulae  were 
derived  included  the  following  functions: 

• Adiabatic  Irreversible  Compression  or  Expansion 

• Isentropic  Acceleration  to  Sonic  Velocity 

• Bivariate  and  Trivariate  Table  Interpolation 

• Mach  Number  and  Pressure  Drop  as  Functions  of  Flow  Function 


3.3  LOW  LEVEL  DERIVATIVE  SET 

In  the  course  of  differentiating  Che  complete  engine  model,  an  alternate 
method  of  applying  the  analytical  derivative  concept  was  devised.  The  devel- 
opment of  this  method  was  carried  forward  on  a parallel  path  under  separate 
funding,  i.e.,  not  Contract  F3361 5-78-C-2203  funds.  As  a matter  of  general 
information  and  to  report  completely  on  analytical  derivative  techniques,  a 
description  of  this  alternate  method  is  included. 
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The  analytical  derivative  technique  for  calculating  the  derivative  data 
needed  to  balance  an  engine  cycle  requires  the  partial  differentiation  of  all 
of  the  component  subroutines.  The  analytical  derivative  method  will  require 
the  differentiation  of  any  new  subroutines  or  as  a minimum,  the  addition  of 
partial  derivatives  to  account  for  any  equation  changes.  Therefore,  time  and 
effort  are  required  to  implement  this  procedure  for  any  new  engine  cycle. 

As  an  alternate  approach,  differentiate  the  lower  level  subroutines, 
such  as  those  subprograms  which  calculate  thermodynamic  and  aerodynamic 
properties,  and  use  these  derivatives  to  update  the  subroutine  output  data. 
This  approach  is  implemented  by  identifying  two  separate  calculation  paths 
in  each  subroutine.  One  of  the  paths  is  the  normal  subroutine  calculation 
and  this  path  is  taken  whenever  a base  point  calculation  is  made.  Added  to 
this  path  are  the  calculations  of  the  partial  derivatives  of  the  output  prop- 
erties with  respect  to  each  of  the  input  parameters,  and  the  storing  of  both 
the  base  point  data  and  the  calculated  partial  derivatives.  This  added  sub- 
program logic  does  increase  the  calculation  time  for  a base  point  (102). 

The  second  path  in  these  subroutines  is  activated  whenever  one  of  the 
independent  variables  is  incremented  for  the  purpose  of  calculating,  by  the 
finite  difference  method,  the  partial  derivatives  used  in  the  iteration 
method  to  balance  the  engine  cycle.  This  alternate  path  calculates  the  dif- 
ference between  the  saved  input  values  and  the  new  input  data  and  then  com- 
bines these  differences  and  the  saved  partial  derivatives  to  calculate  updated 
output  data.  This  path  then  bypasses  the  first  calculation  path  and  returns 
to  the  higher  level  subroutine.  This  second  path  has  considerably  fewer  exe- 
cutable statements  and  therefore  requires  less  computation  time  to  calculate 
an  updated  set  of  output  data. 

This  approach  has  the  advantage  that  once  the  partial  derivatives  have 
been  derived  and  included  in  each  of  the  lower  level  subroutines  they  will 
not  change  and  can  be  used  in  any  computer  cycle  program.  Component  sub- 
routine changes  or  additions  or  a completely  new  component  will  not  require 
changes  in  the  low  level  subprograms  eliminating  the  need  for  updating  the 
derivative  calculation  procedures.  The  time  and  effort  to  incorporate  these 
new  low  level  subroutines  in  any  computer  engine  cycle  program  are  minimal. 

In  order  to  evaluate  this  aproach,  a new  set  of  low  level  subroutines 
was  developed  and  checked  out  for  the  chosen  VCE  cycle.  The  low  level  sub- 
routines which  were  altered  to  include  this  new  procedure  are  listed  below. 

These  subroutines  represent  approximately  102  of  the  total  engine  cycle 
computer  program  logic;  but  since  they  are  called  numerous  times  in  calculat- 
ing a balanced  engine  cycle  operating  point,  they  account  for  about  752  of  the 
processor  time.  Therefore,  any  technique  which  will  lower  the  calculation 
time  of  these  subroutines  will  have  a substantial  effect  on  the  total  com- 
puter processor  time. 
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ADIABX 


I 

Calculates  an  adiabatic  compression/expansion  process. 

Calculates  exhaust  nozzle  throat  area. 

Finds  the  Mach  numbers  and  area  required  by  a variable  area 
mixer  at  the  point  of  confluency. 

Two  dimensional  table  lock  up  procedure  Z » f(X,Y). 

Three  dimensional  table  look  up  procedure  Z * f(X,Y,W). 

Solves  for  Mach  number  given  flow  function. 

Calculates  the  turbine  exit  swirl  angle. 

Calculates  the  temperature  and  the  thermodynamic  properties 
from  enthalpy,  fuel  air  ratio,  and  pressure. 

Calculates  the  temperature  and  the  thermodynamic  properties 
from  entropy,  fuel  air  ratio,  and  pressure. 

Calculates  the  enthalpy  and  the  thermodynamic  properties 
from  temperature,  fuel  air  ratio,  and  pressure. 

Establishes  the  coefficients  required  to  calculate  thermo- 
dynamic properties. 

Calculates  the  dissociated  thermodynamic  properties. 

The  subroutine  MODELX,  which  controls  the  order  in  which  the  components 
are  called,  had  to  be  changed  in  order  to  supply  an  indicator  which  informs 
the  low  level  subroutine  when  a base  point  or  derivative  point  was  being  cal- 
culated. After  incorporating  these  new  subroutines  into  the  engine  cycle 
computer  program  the  calculated  performance  did  not  change  (within  engineer- 
ing significance). 

In  order  to  demonstrate  this  alternate  approach,  one  of  the  low  level 
subroutines  from  the  Air  Force's  computer  program  SMOTE  was  chosen  to  be  dif- 
ferentiated and  reprogrammed.  The  subroutine  chosen  was  PROCOM  because  it 
is  called  a large  number  of  times;  and  if  these  changes  were  implemented,  there 
could  be  a measurable  improvement  in  calculation  time.  PROCOM  calculates 
the  speed  of  sound  (CSEX),  ratio  of  specific-heats  (AKEX) , gas  constant  (REX), 
molecular  weight  (AMW),  specific-heat  at  constant  pressure  (CPEX),  enthalpy 
(HEX),  and  entropy  (PHI),  as  a function  of  temperature  (TEX)  and  fuel-air 
ratio  (FARX).  The  procedure  was  to  differentiate  CPEX,  HEX,  PHI,  and  AMW 
with  respect  to  TEX  and  FARX;  for  example,  for  CPEX. 

CPEX  - (CPA  ♦ FARX  * CPF)/(1.  ♦ FARX) 

Since  CPA  and  CPF  are  functions  of  temperature  only,  then  the  following  is 
the  derivation  of  the  partial  differentiation  of  CPEX  with  respect  to  TEX  and 
FARX. 


L 


AENOZX 

COFFEX 

CURVEX 

CURVSX 

PRXXXX 

SWIRLX 

THOFHX 

THOFSX 

THOFTX 

GETCXX 

HSCALX- 

CSALPX 
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CPA  = + b^T  ♦ c^T^  + d^T^  + e^T^*  ♦ f^T-*  ♦ g^T^  •♦■  h^T^ 


CPF  = 32  + b2T  + c 2T^+  d2T^+  e2^  4+  ^2^  82^  h2T  2 


= CPAUT  = b + 2c  T + 3d  T2  + 4e.T3  ♦ 5f,T4  ♦ 6g,T5  + 7h.T6 
o LEX  ill  l l l l 


1111  = CPFUT  = b2  + 2c2T  + 3d2T2  ♦ ^T3  ♦ 5f2T4  ♦ 6g2T5  ♦ 7h2T6 


3 CPA  = 3 CPF 


3 FARX  3 FAUX 


3CPEX 
3 TEX 


3CPEX 
3 TEX 


3CPEX 


3CpA  . * 3 CPF 

+ FARX  TfEX 


3 TEX 


;]  A- 


♦ FARX) 


= CPEXUT  = (CPAUT  + FARX  * CPFUT)/(1.  + FARX) 


CPF 


CPA  + FARX  * CPF 


3 FARX  (1.  + FARX)  ' (1.  + FARX)2 


substitute  CPA  = (1.  + FARX)  * CPEX  - FARX  * CPF 


3CPEX 
3 FARX 


CPEXUF  = (CPF  - CPEX) /( 1 . + FARX) 


Differentiate  HEX  with  respect  to  TEX  and  FARX. 

HEX  - (HEA  ♦ FARX  * HEF)/(1.  + FARX) 

By  definition 
3 HEX 


3 TEX 


HEXUT  - CPEX 
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and  by  a derivation  similar  to  the  above  for  CPEX 


3 HEX 
3FARX 


HEXUF  » (HEF-HEX) /( 1 . + FARX) 


Differentiate  PHI  with  respect  to  TEX  and  FARX. 

PHI  » (SEA  + FARX  * SEF)/(1.  + FARX) 

By  def init ion 
3 PHI 


3 TEX 

3 PH  I 
3 FARX 


= PHIUT  = CPEX/TEX 


= PHIUF  = ( SEF  - PHI ) / ( 1 . + FARX) 


Differentiate  AMW  with  respect  to  TEX  and  FARX 
AMW  = 28.97  - .946186  * FARX 
3 AMW 


3 TEX 

3 AMW 
3 FARX 


= 0 


-0.946186 


The  remaining  variables  REX,  AKEX,  and  CSEX  are  functions  of  variables 
which  have  been  calculated.  Differentiating  these  variables  and  combining 
base  values  with  partial  derivatives  times  the  change  in  the  independent  var- 
iables leads  to  a more  comlex  expression.  As  an  example,  differentiate  REX 
with  respect  to  TEX  and  FARX. 

REX  = 1 . 986375/AMW  « 1.986375/(28.97  - .946185  * FARX) 

3 REX 


3 TEX 

3REX 
3 FARX 

3 REX  -1.986375  dAMW  1.986375  * 0.946186 


/dREx\  /dAMW  \ 
l dAMW J l^dFARxj 


3FARX  (AMW)2  dFARX 


(AMW) ‘ 


Then 


REX  = [REX] base  + 1 • 879480* (FARX  - [ FARX] base)/ (AMW)  2 
which  is  more  complex  than 
REX  - 1. 986375/AMW 
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In  order  to  calculate  the  updated  information  it  is  necessary  to  save 
the  following  data  when  a base  point  is  being  calculated. 

HIST(INH)  - FARX 

HIST( INH+1 ) - TEX 

HIST( INH+2)  - CPEX 

HIST(INH+3)  = PHI 

HIST( INH+4)  - HEX 

HIST( INH+5)  = AMW 

HIST( INH+6)  - CPEXUT 

HIST( INH+7)  - CPEXUF 

HIST( INH+8)  - HEXUF 

HIST( INH+9)  - PHIUF 

HIST(  INH+10)  =•  LOOPER 

The  updated  data  are  then  calculated  by  the  following  equations: 

DELF  - FARX  - [FARX]baae 
DELT  * TEX  - [TEX]base 

CPEX  » [CPEX]baae  + CPEXUT  * DELT  ♦ CPEXUF  * DELF 

HEX  - i”LX]baae  + [CPEX]base  * DELT  ♦ HEXUF  * DELF 

PLI  * I PHI]base  + [ CPEX]base  * DELT/ [ TEX] baae  «•  PHIUF  * DELF 

AMW  - [AMW] baae  - 0.946186  * DELF 


It  should  be  noted  that  setting  CPEXUT  and  CPEXUF  to  zero  results  in  suf- 
ficient accuracy  to  efficiently  balance  an  AEG  cycle  calculation.  There  are 
two  listings  of  PROCOM  showing  the  current  version  (Table  1)  and  the  version 
with  the  analytical  derivative  logic  included  (Table  2).  In  the  analytical 
derivative  version,  logic  is  included  which  will  assign  the  correct  storage 
space  to  each  call  to  PROCOM  and  which  will  choose  the  proper  path  depending 
upon  whether  the  calculation  is  a base  or  a derivative  point.  Logic  must  be 
added  to  the  ENGBAL  program  to  set  the  indicator  IBASE  to  the  correct  value 
depending  on  whether  a base  or  a derivative  point  is  being  calculated.  The 
analytical  derivative  version  assumes  that  the  variable  LOOPER  is  increased 
on  every  path  through  the  engine  cycle  calculation  and  that  IBASE  will  be  set 
equal  to  LOOPER  whenever  a base  point  calculation  is  made.  In  addition,  the 
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Table  1,  Listing  of  Subroutine  PROCOM  (Thermodynamic  Property  Calculation). 


1190 

X 

1200 

X 

1210 



1220 

X 

12  30 

X 

1240 

X 

12  50 

1260C 

12  70 

12  80 

X 

12  90 

X 

1 300 

1310 

X 

1320 

X 

1330  - 

_x 

13  40 

13  50 

X 

1360 

X 

SUBROUT 1NE  PROCOM(FARX*TEX#CSEX.AICEX,CPEX,fiEX#PHI#HEX) 

If  (FARX.LE. 0.067623)  60  TO  1 
FARX-0. 067623 

WR  I TE ( 8# 1 01  ) . 

1 If  (TEX.  6E.  300.  1 6 C TO  2 
TEX-300. 

WRITE(8#102>  . 

2 If  (TEX. IE.  4000.)  60  TO  3 
TEX-4000. 

WRITE  (8*103)  _ ... 

3 If  (fARX.GE.0.0)  GC  YO  4 
FARX-0.0 

_WR  1TE  < 8#  1 04J  _ _ . 

AIR  PATH 

4 CPA  «({ ( ( (<1  ,0115540E-25*TEX-1  .4S26770E-21  1 * T EX 

X *7.6215767E-1  81  * TEX-1 . 51 2 82 59E-1  4>*TEXr6.7J  7 8376E-1  21  * T EX 
X *6.551 9486E-081*TEX-5.1 5 36879E-0 5 > » T EX* 2 . 50 20051 E-01 
HEA*(C( ( ( ( (1 .2644425E-26*TEX-2. 0752522 E-22)*T EX 
♦1.2702630E-1 81* TEX-3. 025 651  8E- 151  * T£X-1. 6Z 94594 E-1 21 « I E X 
*2.1 83 982 6E-081* TEX-2. 5768440E-0 5 ) * T EX*2 . 50 2 00 5 1 E-01 > * TE X 
-1.7558886E*00 

SEA»*2. 50  200  51  E-01*  ALOG<  TEX1  ♦((<(<  41 . 4450  767  E=l26*T  EX 
-2.4211 288E-22)*TEX*1 .5243153E-1 8)* TEX-3. 7820648E-1 5)* TEX 
-2.2  392  7 90E-1 2 1 * TEX*3. 2 7 5 9 743E-081 * T E X-5 . 1 5 7 6879E -051  *TEX 

♦ 4.5  43-2  3QQEr  0 2 . 

If  (FARX.LE. 0.01  60  TO  5 
FUEL/ AIR  PATH 

CPF  «<<<4.4{?.2678710E-23*TEA-J.333566AE=201*X£X 

♦1.021 291 3E-161*TEX-4.2051104E-131*TEX*9.9686793E-101*TEX 
-1.3771901 E-0 63*TEX*1. 22586 30 E-03)*TEX*7.3816638E-02 

HE  F«XU.a.(  (9^084  8388£T26*-IEAr.U.9II5fljL4  9£t^UaIXX 

♦1.7021 52 5E-1 71* TEX- 8.41 02208E-1 41 *TEX*2. 4921 698E-101* TEX 
-4.5906332E-071*TEX*6.1293150E-04)*T£X*7.381 6638E-02) 

.*TEX*3.0  581530£t01 ... 

SEF»*7.381 6638E-02* AL  06(TEX)*((((((1.0382670E-25*TEX 
-2. 222611 8E -2 1)*TEX*2. 0425 82 6E-17)*TEX-1, 051277 6E-13)*TEX 


1370  X *6.4 83398E-01 

1380  5 CPEX*<CPA*FARX*CPFJ/<1.*FARX1 

-1310 MiX*  l HE  A*  F ARX  « H£f l/(  1 . ♦ FARXl 

1400  PHI-(SEA*FARX*SEF1/<1  . *FARX) 

1410  AHW-28.97-.9461 86*FARX 

_1420 REX-J  .986375/JULk 


1430  AKEX-CPEX/ (C  PEX— RCX1 

1440  CSEX-SQRT(AKEX*REX*TEX*2S031.37) 

.1430 REJUR  N 


1460  101  FORMAT(1h0#63HINPUT  FUEL-AIR  RATIO  ABOVE  LIMITS#  IT  HAS  BEEN  RESET 

1470  X TO  0.067623-6HSSSSSS) 

1480 102 — F0AHAlIXH0#35HgJ0CIlJI.-lNPUT  IENPE-BAIU1E  BELOW  3Q0w6H**AAAS1 


1490  103  FORHAT(1hO#3AHPROCOM  INPUT 
1500  104  EORMAT(1hO#38HPROCOM  INPUT 
1510 EBB 


TEMPERATURE  ABOVE  4000. #6HSS*SS*> 
FUEL-AIR  RATIO  BELOW  EER0#6MSSSSSS> 
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-Table  2.  Sample  Listing  of  Subroutine  PROCOM  (with  Analytic  Derivatives) 
(Subroutine  Not  Checked  Out). 


1000 SUBROUTINE  PROCOM(FARXsTEX«CSEXsAKEXs.  CPE-X/REXsPHlsHEX.lNH) 

1010  C0MM0N/HISTRG/I8ASE,INHS#L00PERsHIST(891  ) , INHK81  ) 

1020  DIMENSION  IHIST(891> 

.1030  ..  _ EQUIVALENCE  ( H I S T <11#  I H LS  T 41  J ) 

1040  IF  (INH.NE.0)  GO  TO  10 

1050  I N H* I NH  S 

J060 INHS*INHS*11  ... 

1070  10  IF  (IHIST ( INh  + 10).NE.IBASE)  GO  TO  20 

1080  DELF»FARX-HIST<INH) 

1090.  ..  DELT«TEX-HIST(INH*1)  . 

1100  CPEX»HIST( INH*2  >*HIST (INH*6)*DELT*HI ST(INH*7)*DELF 

1110  HE  X*H  I S T { INH+4 ) *HIST<  I NH ♦ 2 ) * DE L T*H I S T ( INH ♦ 8 ) * DEL F 

1120 PHI«HISTIINH+3lTHlSTCINH+2J*D£i.T7HISJXUmAU ... 

1130  X *H1ST( INH*9)*DELF 

1140  AMU*HIST(INH*5>-0l946186*0ELF 

1150 GO  TO  6 

1160  20  IF(FARX.LE. 0.067623)  GO  TO  1 

1170  FARX*0. 067623 

11.80 —WRITE!  SHOD 

1190  1 I F (TEX. GE. 300. ) GO  TO  2 

1200  TEX*300 . 

1210 -WRITE  (8*102)  — 


1220  2 IF(TEX.LE.4000.)  GO  TO  3 

1230  TEX»4000. 

1.2 AQ WILII2{8.,103X 

1250  3 IF(FARX.GE.O.O)  60  TO  A 

1260  f A RX*0 . 0 

1220 U R XTE  48  .1  0 41 

1280C  AIR  PATH 

1 290  4 CPA  ■ ({((  ((1  .01  1 5 54 0E-25*TEX-1  . 4526770 E- 21  )*TEX 

1300 X_*-7.621  57  67E.-1 81  ■>  J2X=3.5323259£-t51-»3EX-6.  71  7t37-6£=12l-«IEX 

1310  X ♦6.551 94866-08>*TEX-S.1 536879E-05) * T EX*2. 5020051 E-01 

1320  CPAUT  »((<((7. *1.0115 540E— 2 5*TEX-6.* 1 • 4526770E-2 1 ) *TEX 

1 330 X~*5.*-7.4213767E-18)*TEX-4.*.U512S259E=14TL*72-X=2.Jt4.71783  76E-12)  - 

1340  X *TEX*2.*6.5519486E-08)*TEX-5.1536879E-05 

1350  HEA*( (( ( (( (1.2644425E-26*TEX-2.0752522E-22)*TEX 

1380 -X— * 1.2  7026306^1  81*  TEX-3^0 256318 E--1 51*  TEX.=-1  .67945942=1 21 « TEX 

1370  X *2.1 839826E-08)*TEX-2. 5768440E-0 5 >* T EX*2. 5020051 E-01 > *T EX 

1380  X -1 .75588866*00 

1390 $£A«*2.5020051£'=01*AL062T-2Xl-*.X4.<-<.2tl-.445Q7-622=,25*-I2X 

1400  X -2.421  1 2 88E-22)*TEX*1. 52431 53E-1 8)* TEX-3. 7820648E-1 5) *TEX 
1410  X -2.2 392 7 90E-1 21* T6X*3. 2 759743E-08)* TEX-5. 15 788796-051* TEX 

1420 x_* 4.54323002=02 

1430  !F<fARX.LE.O.O>  GO  TO  5 

1440C  FUEL/AIR  PATH 

1450 CP-f-J*(t  (((X7.267871 02=23 *X£X-1 . 33356686-20  )*7£« 

1460  X *1.021  291  3E-1  6)*  TEX-4.  2051 104  E-1  3)  *TEK*9.9686793£-10W£X 
1470  X -1.3771901 E-06)*TEX*1 .2258630E~03>  * TEX* 7. 381 663 8 £-02 

1480 C-P-FUJ  «22U  (7.-*7.  26757-1 0£»25*J£X«6.*  1.353  58888*20 J *TEJC 

1490  X *5. *1.021291 3£-1 6) * TEX-4. *4. 2051 104 £—13) *T£X*3« *9.96867 93E-10) 
1500  X *TEX-2.*1.3771901£-06)*T£X*1. 22586306-03 

1520  X *1.7021 525E-17J*TEX-8.4102208£-14>*TEX*2,4921698E-10)*T£X 

1530  X -4. 5 906 3 32E -0 7 >*T£X*6. 12 931 506-04) *TEX*7 .38166386-021 

1540 X • I2X*3.Q 5815302*01 

1550  SE  F»*7. 381 6638E— 02*AL06<TEX)*< (< ( <<1 .03826706— 25*T£* 

1560  X -2.22261 1 BE-21 )*TEX*2.0425826E-1 7)*T6X-1 .05 1 2 7766-1 3) *T EX 
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THIS  PAGE  IS  BEST  QUALITY  PRACTICABLE 
PROM  COPY  FURNISHED  TO  DOC  - — ' 


Table  2.  Sample  Listing  of  Subroutine  PROCOK  (with  Analytic  Derivatives) 
(Subroutine  Not  Checked  Out) (Concluded) . 


— 1570  . X ♦ 3.3228928E-10)*TEX-6.88S9505E-07)«TEXAl..2238630E-03J*rEX 

1580  X ♦6.483398E-01 

1 590  5 DUM«1 . 0 ♦ F ARX 

..  1600  CP£X«UPA*FARX*CPf)yDUM  . . 

1610  CPEXUT»(CPAUTeFARX*CPFUT)/DUM 

1620  CPEXUF»<CPF-CP£X)/DUM 

1630  _ . -HE.X»IHEA*F  ARX*HEF)7DUM 

1640  HEXUF*(HEF-HEX)/DUM 

1650C  MEXL’T«CPEX 

_ 1660  ...  PHI«(SEA*FARX*S£F)/BUM  

1670  PHIUF»<SEF-PHI)/DUM 

1 6 80  C PHIUT*CPEX/TEX 

16  90  AMM*28..97-.9461  86*/ARX 

1700  HI ST< INM) »FARX 

1710  HI ST( INM*1 )»TEX 

1720 HI ST( 1 NH*2 ) *CPEX 

1730  HISTt INH»3 )«PHI 

•740  HI $T(INH»4)»H£X 

1730  HIST.11NH45)»AHW . 

1760  HI$T< INH*6)»CPEXUT 

1770  HIST( INH*7>«CPEXUF 


1780 HIST(1NH*8)«HEXUF 

1790  MI ST< INH+9)«PHIUF 

1800  lHlST<lNHdO)«LOOPER 

1810 6 RE  X»1_. .986.3 75  / AMW . 

1820  AKEX«CPEX/(CPEX— REX) 

1830  CSEX«SaRT (A<EX*REX*TEX*25031 .37) 

1840 RETURN _ _ 

1 850  101  F0RMAT(1  H0.63HINPUT  FUEL-AIR  RATIO  AB0\/E  LIMITS.  IT  HAS  «EEN  RESET 

1860  X TO  0. 067623. 6HSSSSSS) 

1870 10 2 — f .OiLM ALLLH Q-.J. 5.MPR0C0M_JLN& U J_J_£HP £ RA.T UR.£_BJLLOM 300 . . 6HSSSSSS)  

1880  103  FORMAT (1  HO. 36HPR0C0M  INPUT  TEMPERATURE  ABOVE  4000. .6 HSSSSSS) 

1890  104  FORMATd H0.38HPR0C0M  INPUT  FUEL-AIR  RATIO  BELOW  2ER0.6M* SSSSS) 

...  1.90.0 E N P 


call  sequence  Co  PROCOM  must  be  increased  Co  contain  a unique  history  array 
pointer.  Since  there  are  81  possible  calls  to  PROCOM,  the  history  storage 
array  (HIST)  is  dimensioned  by  891  (81*11)  and  the  history  pointer  array 
( INHI)  by  81. 

There  is  added  logic  at  the  beginning  of  the  subroutine  PROCOM  which 
determines  if  that  specific  call  to  PROCOM  had  been  executed  previously  and 
if  not,  sets  the  storage  indicator  (INH)  and  updates  the  history  counter 
INHS  by  11.  This  logic  has  been  added  to  ensure  that  the  first  time  that 
each  call  to  PROCOM  is  executed  the  history  storage  array  for  that  storage 
indicator  has  been  filled  with  correct  data. 

This  alternate  approach  was  incorporated  into  the  VCE  cycle  and  the 
performance  for  the  411  approved  flight  operating  points  was  calculated. 

The  results  showed  that  this  approach  used  47%  less  computer  processor 
time  than  the  original  computer  program  but  required  5K  more  computer  stor- 
age. The  alternate  method  is  easily  implemented  using  minimal  time  and  cost. 


4.0  APPLICATION  OF  THE  METHOD  TO  SMOTE 


4.1  CONTROL  LOGIC 

As  a vehicle  to  illustrate  the  application  of  the  method  using  the  com- 
plete analytic  derivative  set,  the  Air  Force  engine  cycle  simulation  program 
was  selected.  This  program,  titled  SMOTE  (Simulation  of  Turbo  fan  Engine)  was 
developed  in  the  Components  Branch,  Turbine  Engine  Division,  Air  Force  Aero 
Propulsion  Laboratory  and  is  described  in  References  2 and  3.  A simultaneous 
Newton-Raphson  iteration  method  is  used  to  calculate  balanced  cycle  perfor- 
mance. The  partial  derivatives  are  obtained  by  numerical  differentiation. 

A matrix  of  differential  error  equations  is  then  solved  to  determine  the  cor- 
rect values  of  the  independent  variables  which  would  produce  zero  errors . A 
flow  chart  of  the  program  is  shown  in  Figure  2. 

As  an  example  of  the  method,  the  derivation  of  the  analytical  derivatives 
and  the  associated  computer  program  logic  for  a compressor  and  a combustor 
will  be  shown.  It  should  be  pointed  out  that  none  of  the  Fortran  IV  sub- 
routine listings  referred  to  in  this  section  have  been  checked  out.  The 
listings  are  intended  as  an  aid  in  understanding  the  examples,  not  as  execu- 
table code.  Control  of  the  SMOTE  program  ia  contained  in  the  main  subroutine 
ENGBAL.  This  subroutine  controls  all  engine  balancing  loops,  checks  toler- 
ances and  number  of  loops,  and  loads  the  matrix.  The  functioning  of  this 
subroutine  is  shorn  schematically  in  Figure  2,  and  a listing  of  the  sub- 
routine is  given  in  Table  3.  The  incorporation  of  analytic  derivatives 
would  necessitate  the  rewriting  of  the  ENGBAL  subroutine.  The  new  flowpath 
would  be  similar  to  that  shown  in  Figure  3.  In  Figure  3,  the  differentials 
of  the  independent  variables  are  designated  by  the  symbol  VU  with  an  appro- 
priate subscript.  These  differentials  are  incremented  sequentially  (from  1 
to  6)  by  unity.  When  the  indicator  IDERV  ■ 1,  the  deck  is  on  a derivative 
path.  In  this  case,  the  subroutine  base  point  logic  (the  block  labeled 
ENGINE)  is  jumped  and  only  the  derivatives  are  calculated  in  the  individual 
subroutines.  The  errors  can  be  calculated  from  the  equation: 

ERR(I)  - ERR(  I)\  ^SE  ♦ d (ERR(I)) 

7 POINT 

The  error  differentials  are  obtained  by  differentiating  the  individual 
error  definitions.  The  unknown  terms  will  then  represent  derivatives  and/or 
differentials  evaluated  in  the  individual  subroutine. 

For  example.  Table  4 shows  the  independent  and  dependent  variables 
(errors)  for  the  SMOTE  deck.  The  differentiation  of  the  first  error  would 
proceed  as  follows: 
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— k/1, 

Table  3.  Listing  of  Subroutine  ENGBAL. 


QUALm  PR4CI1CABLI 
raaM  copX  BURNISHES)  TO  DOC 


1000  SURROUTINF  FNGRAL 

10 1 0 COMMON  / ALL/ 


1 0?0 

1 WORD 

, IDF5  , JDFS 

.KPFS 

, MODF  , I NIT 

. IDUMP  . iaktp  . 

1030 

2IGASMX 

, IDBURN, IAFTBN 

. IDCD 

, IMCD  , IDSHOC 

, IMSHOC.N07FLT. 

1040 

3 1 TRY  S 

,t.OOPFR,  NOMAP 

.NUMMAP 

, MAPFDG.TOI  ALL 

,FRP( 4) 

1050 

COMMON 

/DESIGN/ 

1060 

1PCNFGU 

.PCNCGU.T4GU 

. DUMD  1 

, DUMD?  .DELFG 

, PF1.FN  .DELSFC, 

1070 

2ZFDS 

.PCNFDS, PRFDS 

. FTAFDS 

.wafds  .prfcf 

.ftafcf.wafcf  . 

1080 

3ZCDS 

.PCNCDS, PRCDS 

. ETACDS 

.WACDS  .PRCCF 

. FTACCF, WACCF  , 

10Q0 

4T4DS 

,WFPDS  .DTCODS 

.FTARDS 

.wascps.dpcops 

, TTCOCF , FTA  RC  F , 

1 1 00 

5TFHPDS 

.CNHPDS, ETHPDS 

, TFHPCF 

.CNHPCF.ETHPCF 

,DHHPCF,T2DS  . 

1 11  0 

4TFLPDS 

.CNLPDS.FTLPDS 

.TFLPCF 

.CNLPCF, FTLPCF 

, DHF..PCF.T71  DS  . 

1 1 ?0 

7T24DS 

, WFDDS  , DTDUDS 

.ETADDS 

. WA  23DS , DPDU  DS 

, DTDUCF , ETA  DCF , 

1130 

8T7DS 

, HFADS  .DTAFDS 

.ETAADS 

,.WG4CDS,DPAFDS 

, DT  AFCF , FT  A ACF , 

1 1 40 

9A55 

,A26  , A6 

,A7 

,A8  ,A9 

.A2R  . A ?9 

1 1 50 

APS55 

,AM55  ,CVDNOZ 

, CVMNOZ 

, A8SAV  , A9SA  V 

, A28SAV.A29SAV 

1140 

COMMON 

/ FRONT/ 

1170 

IT  1 

,P1  ,H1 

.SI 

,T2  ,P  2 

,H2  ,S2 

1 1 R0 

7T7 1 

.P21  , H2 1 

. S2 1 

,T3  ,P3 

,F3  .S3 

1 190 

3T4 

.P4  ,H4 

,S4 

«T5  . ,P5 

,H5  ,S5 

1200 

4T55 

,P55  ,H55 

, S55 

, BLF  , BLC 

, PF.PU  ,RLOR  . 

1210 

5CNF 

, PRF  , ETAF 

. WAFC 

,WAF  , WA3 

,h'G4  .far  4 , 

1220 

6CNC 

,PRC  f ETAC 

, HAOC 

. NAC  fETAt3 

, DPCOM  .DUMF  , 

1 ?30 

7CNHP 

, FTA  THP . D HTCHP 

,DHTC 

, BLHP  ,NG5 

. FAR5  .CS 

1240 

8CNLP 

.ETATLP.DHTCLP 

,DHTF 

• BLLP  , WG55 

, FAR55  .HPFXT  , 

1?50 

9 AM 

, A LTP  , FTAR 

. ZF 

, PCNF  ,ZC 

. PCNC  ,WFB 

1 260 

ATFFHP 

.TFFLP  . PCBLF 

. PCBLC 

.PCBLDU.PCBLOB 

, PCBLHP , PCBIJ.P 

1770 

COMMON 

/ SIDE/ 

1 280 

XXP1 

, XWAF  , XWAC 

. XBLF 

. XBLDU  , XH3 

.DUMSI  .DUMS2  . 

1290  ~ 

XXT21  ' 

,XP21  " , XH2.1 

. XS2 1 

,T23  _ .P23 

,H?3  ,S73~  , 

1300 

3T24 

,P24  ,H24 

,S24 

,T25  . P25 

. H25  ,S25 

1310 

4T28 

«P28  ,H28 

. S28 

,T29  ,P29 

« H29  ,S29 

1 370 

5WAD  “ 

,WFD  " .WG74 

,FAR?4 

, FTAD  .DPDUC' 

. BYPASS, DUM53  . " 

1330 

6TS2R 

.PS2R  ,V28 

, AM28 

,TS29  ,PS29 

• V29  , AM29 

1340 

COMMON 

/ RACK/ 

'1350 

""  XXT55 

XP55  ,XH55 

. XS55 

« XT  25"  .XP25 

, XH25  ” '.XS75  ". 

1360 

XXNFB 

.XWG55  , XFAR55 

, XWFD 

.XHG74  .XFAR74 

,XXP1  .DUMB  , 

1370 

3T6 

, P6  ,H6 

• S6 

»T7  ,P7 

.H7  ,S7 

1380 

4T8 

,P8  ,H8 

,S8 

.T9  , P9 

'.H9  ““  ,S9  , 

1390 

5WG6 

,WFA  ,WG7 

.FAR7 

. FTAA  .dpaft 

. V55  ,V25 

1400 

6PS6 

,V6  , AM6 

,TS7 

,PS7  ,V7 

, AM7  , AM25  . 

141  0 

7TS8 

,PSR  ,VR 

. AM8 

,TS9  ,P59 

, V$  , AM9  , 

14  20 

8VA 

,FRD  , VJD 

,FGMD 

. VJM  ,FGMM 

, FGPD  ,FGPM  , 

1430  9FGM  ,FGP  ,WFT  ,WGT  .FART  ,FG_  . FN  ^SFC 

1440  “ DIMENSION  VA R< 6 ) ", DEL < 6 ) , ERRB< 6J", DEL  V AR(~6 f ,PMAT ( 6 , 6j , Vm  AT (ft ) . 

1450  I AMAT ( 6 ) 

1440  DATA  AW0RD/4HFNGBAL/ 

14 70  DATA"  VDELTA’,  VlTM • VCHiNCf sx/ 

1480  1 0.001 ,0. 100,0.850,4/ 

1490  CALL  INPUT 

'1500  IFf  fNIf. EO.T)  "GO 'To  50  "" 

1510  TFFHP-TFHPDS 

15  70  TFFLP-TFLPDS 

J530#50  LOOPlRi’0  ' 

1540  NUMMAP-0 

1550  NOMI SS“0 

1540#  1 LOOP-0 
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1*7* 

15  8* 
1590 
1 *** 
1610 
1620 
1630 
1640 
1 4**#? 

I 660#3 
1 *7* 
1680 
1690 
I 700 
1710 
1 720#55 
1730 
1 740 
1750 
1760 
1770 
1780 
1790 
1800 
18  10 
1820 
1830 
1840 
I850#5 
i860 
1870 
1880 
1 800 
190046 
1910 
I9?0#T 
1930 
1940 
1950 
I960 
1970 
I 980470 
1990 
2000 
"2010 
?0?0#8 
203049 
" 2040 
2050 
2060410 

2070 

2080 

2090 

' “21 00  - - 
2110 
2120 
' 21  30" 


Table  3.  Listing  of  Subroutine  ENG UAL  (Continued). 


THIS  PAGE  IS  BEST  QUALITY  PRACTICABLE 
JBOM  COPY  EURW1SHED  TO  DBC  


MISMaT-0 
NOMAP-0 
IGO-2 

no  2 1-1,6 
VMAT(I)-0. 

AMAT ( I )-0, 

HFI.VARd)-*. 

DO  2 L-1,6 
FMAT(  I ,U  -0. 

LCX)PER«LOOPER+ 1 
CALL  COFAN 
WORD-AWORD 

IF(L(X)PER.GT.  ITRYS)  GO  To  20 
IF(NOMAP.GT.0)  GO  TO  1 
NUMMAP-0 
VAR( 1 )*ZF*100. 

IF(  MODE.  NE  .3)  VAR(2)»PCNF 
IF(M0DE.E0.3)  VAR(2)*T4/10. 

VAR( 3)«ZC* 100. 

IF f MODE. NE. 1 ) VAR(4)«PCNC 
IF(MODF.FQ. 1)  VAR(4)«T4/10. 

VAR(5;-TFFHP 

VAR(6)-TFFLP 

DO  4 1-1,6  . - . 

IF(ABS(ERR( I )) .GT.TOLALL)  GO  TO  5 
i CONTINUE 

CALI  PERF - 

CALL  FRROR 

IF(LOOP.CT.0>  GO  TO  7 

MAPEDG-0 ' - 

MAPSET-0 
DO  6 1-1,6 

FRRB ( I )»FRR( I > 

DEL(I)-VDELTA>VAR(I) 

GO  TO  9 

" IF(MISMAT.GT.0J  GO  TO  30  " " " 

IF(MAPEDG.EO.0)  GO  TO  70 
MAPEDG-0 

MAP 'SET*  I ‘ “ 

VAR  ( LOOP ) -VA  R(  LOOP) ♦ 2.  *DEL  ( LOOP) 

GO  TO  10 

' IF(MAPSET.EQ.0)  VAR(E00Pj-VAR(L(X)P)'eDEL(L(XiP) 
IFfMAPSET.  EO.  1 ) VA  R ( LOOP )— VAR ( LOOPJ-DEL  (L(X)PJ 
MAPSET-0 

■ ffinrT«f,v  “ ' 

EMA T(  I , LOOP )- ( E RRB  < I )-E  RR ( I ) ) /DEL  ( LOOP ) 
LOOP-LOOP* 1 

"I FfLOOP.GT. 6)  "Go  TO  IT  " 

VAR  ( LOOP ) - VA  R(  L OOP  > -DEU  LOOP ) 

ZF-VAR1  U/100. 

‘ IFfMODE.  NE  .3 ) "PCNF-VAW  2 J 

IF(MODF. E0.3)  T4-VARX 2) *10. 

ZOVAR(3)/100. 

*"TF(MODE".NE.  1)  PCNC-VABUJ 

IFfMODE.EO.I  ) T-4— VAR(4J*I0, 

TFFHP-VARX5) 

TFFLP-VaR( 6T 
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Table  3.  Listing  of  Subroutine  ENG BA ~ (Continued). 


2140  IF(ZF.LT.O.)  ZF-0.05 

2150  IF (ZC. LT.P . ) ZC-0.06 

2I60  GO  TO  ( I ,3), ICO 

21  7(1#  I I DO  I 2 1*1,6 

2 I 8 0#  1 2 AMATU  )*-ERRB(I  ) 

219(1  DO  1 d 1*1,  A 

22 OO  IZERo-0 

2210  no  13  LOOP* 1,6 

2220# I 3 IF(PMAT<I.LOOP).E0.0.)  I ZFRO-I ZERO# I 

2230  IF( IZEPO.I  T.6)  GOTO  |4 

2940  HRITM  A,  100)  I 

2250  L(X)PER*ITRYS#  100 

22 AO  GO  TO  20 

2270# 1 4 CONTINUE 

2280  DO  16  LOOP* i, 6 

2290  IZFRO*0 

2300  DO  15  1*1,6 

23 1 o#i5  IF(EMaT(I  .LOOP)  .FQ.O.)  I Z FRO* I ZERO# I 
2320  IF( IZERO.LT. 6)  Go  To  16 

2330  NRITE(6,  101)  LOOP 

2340  LOOPFR=ITRYS#1O0  ~ ' 

2350  GO  TO  20 

9 21A0#1A  CONTINUE 

2370# 1 7 CALL  MATRIX (EMAT, VM AT, AM AT) 

9380  LBIG*0 

2390  VARBIG-0.0 

2400  * DO  18  L*  1,6 

2410  *BSVAR«ABS(VMAT(L)) 

2420  IF( ABSVAR.LE . VLI W*VAR(L) ) GO  TO  18  _ 

2430  ' IF(ABSVAR.LF.VARRIG)  GOTO  18  ' ' 

2440  I.BIG-L 

2450  VARBIG-ABSVAR 

24 AO# 1 R CONTINUE 

2470  VRATIO-1.0 

24R0  IF(LBIG.GT.O)  VRATlO»VLIM*VAR<LBIG)/VARBIG 

24  90  ERRAVE-0.0 

25e0  VMTAVE-0.0 

2510  _ DFLA VF*0. 0_ 

2520  DO  ' 19  L*  1,6  -----  - 

2530  0ELVAR(L)*VRATIO*VMATtL) 

2540  FRRAVE«ERRAVE#ABS<AMAT(L>) 

2550  VMTAVE*VMTAVE#ABS( VMAT(L) ) ' ' “ 

2560  DELAVE-DELAVE#ABS(DELVAR(L)) 

2570# 19  VAR(L)*yARl L) ♦DELVARfL) 

25R0"  "FfiRAVF-ERRAVE/6. " 

2590  VMTAVE-VMTAVE/6. 

2A00  nFLAVF-nPLAVF/6. 

2610  I F ( MTSMAT. ST.OT  Go  To"3 2 

26  20  IF ( NOMI SS . EO .*0)  MISMAT-1 

2630  IFtMISMAT.EO.0J  IGO- 1 

‘ 2640*20  " WR7  TEf  871 02  riOOPER  “ 

2650  DO  21  I*  I, A 

2660#2 1 HRITE(8, 103)  AMAT( I ) , ( EMAT( I ,L ) ,L* 1 ,6) , VMAT( I ) , DELVA R( I ) ,VAR( I ) 

‘ “ 2A70  HRTTF/R",  104)'FRRaVE,VMTAVF.DFLAVF  * • 

26R0#22  IF( LOOPER.LT. ITRYS)  Go  TO  10 

2690  CALL  ERROR 

"2700  'RETURN 
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27I0#30 
2720 
2730#3  I 
2740 
?7R0#32 
2760 
2110 
2780 
2790 
2800 
2RI0 
2820 
2R30 
2840# I 00 
2880# 1 0 1 
2860# I 02 
2970 
2880# I 07 
2890# I 04 
?O00# | 08 
2910 
29  20 
2930 

2940# |06 
29  90 
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Table  3.  Listing  of  Subroutine  ENGBAL  (Concluded). 


VMTAVX*VMTAVE 
DO  3 ! 1-1,6 
AM AT ( I ) =-FRR(  I ) 

GO  TO  17 

WRITRfl,  H*)  AMAT.nRRAVC.nFI  VAR.nFTAVP.VMAT.VMTAVF.VAR 
M ISMAT*M ISMAT+ I 

IF( VMTAVE.LT. VCHNGE*VMTAVX)  GO  TO  2? 

WRITFCR. 106) 

IF<  MISMAT.LT. NOM I SX)  NOMISS-I 

MISMAT=0 

I.(X)P*0 

I GO*  2 

GO  TO  89 

FORMA T(4H0ROW,  12.  I6H  IS  ZERO  IN  EMAT) 

FORMAT!  7H0COLUMN.I2,  16H  IS  ZFRO  IN  FMAT>" 

FORVATC8HB  ERRB.28X23HERROR  MATRIX  AFTER  LOOP  , 14 , 29X4HVMA 
16X6HDELVAR.7X I4HVARI  ABLESSSSSS) 

FORMAT( IH0.FR.4.RX6FI0.4, 10XF10.4.F 1 1.  4.4XF I 1 .4.6HSSSSSS) 
FORMAT(  IH0,F8.4,32X 14HAVERAGE  VALUES, 3 I X, 2F  1 1 .4 , 6HS SS SSS ) 

FORMAT! 1 2H0 AMAT, 7F 16. 6, 6HS SSSSS , _ 

I/,  12H -DELVAR,7F 1 6. 6,6HS$SS SS, 

2/,  1 2H VM AT, 7F 1 6. 6, 6HSS  SS  SS. 

3/,  1 2H VAR,6F!6.6,6HSSSSSS) 

FORMAT!  1H0, 50X 22HCHANGE  TOO  SMAIJLSSSSSS ) 

FND 


Figure  3.  SMOTE  Computer  Program  Flow  Chart  (with  Analytical  Derivatives). 


TFFLP 


ERR(l)  - (TFHCAL  - TFFHP)/TFHCAL 

dERR( 1 ) - TFHCAL  * (TFHCAL  - dTFFHP)  - (TFHCAL  - TFFHP)  » d TFHCAL 

(TFHCAL)2 


dERR(l) 


dTFFHP  + TFFHP  ♦ d TFHCAL 
TFHCAL  (TFHCAL)2 


dERR(l)  - JlIBLS  */dTFFHP\  /TFFHP  V ZdTFHCAlN 
VTFHCAl/  \TFFHP  y \TFHCAL j V TFHCAL J 


dERR(l)  - + 


/ TFFHP  \ * [7dTFHCAL\  /dTFFHPN 
\TFHCA LJ  [^TFHCAL  ) ~ yTFFHP  ) 


Note  that  all  terms  on  the  right-hand  side  of  the  above  equation  would 
be  evaluated  at  the  base  point.  The  quantities  TFFHP  and  TFHCAL  are  normally 
computed  when  the  base  point  is  calculated.  The  two  differential  terms  are 
evaluated  on  the  derivative  path. 

In  general, 

dERR(l)  - -ET^(1)dV.  ♦ } dV  ♦ ...+  dv 

3Vj  1 3V2  2 3Vn  n 

On  any  given  derivative  path  all  the  dV's  will  be  zero  except  one  whose 
value  will  be  unity. 

Therefore , 

dERR(l)  * ^ ^ after  the  first  pass  (L  ■ 1) 

dERR(l)  ■ after  the  second  pass  (L  ■ 2) 

and  dERR(n)  ■ after  the  nth  (last)  pass. 

dVn 


The  subroutine  ENGBAL  as  it  might  appear  with  the  analytic  derivative 
control  logic  in  place  is  shown  in  Table  5.  The  subroutine  AERR  referred  to 
in  the  listing  would  return  the  individual  error  derivatives,  i.e.,  dERR(i)  - 
ERRU(I).  It  would  contain  the  logic  for  all  the  error  term  derivatives  de- 
rived in  a manner  similar  to  that  shown  for  the  first  error  term.  The  base 
point  values  and  the  value  of  the  derivative  and/or  differential  terms  could 
be  communicated  to  the  subroutine  via  labeled  common. 
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Table  5.  Sample  Listing  of  Subroutine  ENGBAL  (with  Analytical 
Derivative  Logic)  (Subroutine  Not  Checked  Out). 

1000  SUBROUTINE  ENGBAL 

10  10  COMMON  / All/ 

’*->0  IWOHD  , I DFS  .JOES  .KPFS  .MODE  .INIT  . IDUMP  , I AMTP 

1030  PIGASMX.IDBURN.IAFTBN.IDCD  ,IMCD  , IDSHOC, IMSHOC.NOZFLT 

I0-J0  3 1 TRY  S ,!.(X)PrR, NOMAP  .NUMMA  P,  MAPFDC  ,TOI  A LI  , FRR(  A) 

1050  COMMON  /DESIGN/ 

1/60  IPCNFGU.PCNCGU.T4GU  .DUMDl  ,DUMD7  .DELFG  , DEI.FN  .DELSEC 

1070  7ZFDS  , PCNFDS, PRFDS  . FTAFDS. WAFDS  .PRFCF  . FTAFCF . WAFCF 

1080  3ZCDS  ,PCNCDS,PRCDS  . ETACDS, NACDS  ,PRCCF  , ETACCF. NACCF 

1->C30  4T4DS  , WFBDS  , DTCODS,  ETaRDS,  WA3CDS,  PPCoOS.rTCOCF  , FTAPCF 
U00  5TFHP  DS.CNHPDS.ETHPDS.TFHPCF, CNHPCF .ETEPCE,  DHHPCF,T2D5 

M > 0 aTFI.PDS.CNI  PDS  , ETLPnS.TFLPCF.CNI  PCF, FTIPCF,  PFLPCF ,T2!  DS 

1120  7T24DS  , WFDDS  , DTDUDS. ETADDS, WA 23DS. DPDU DS, DTDUCF, ETA  DCF 


1 ! 30 

8T7DS 

, WFADS 

.DTAFDS 

.ETAADS 

, WG6CDS 

.DPArDS 

.PTaFCF.FTaaCF, 

1 1 40 

9 A 55 

, A 25 

, A6 

.A7 

,A8 

. A9 

.A  28  .A ?9 

1 1 50 

APS55 

, AM55 

.CVDNOZ 

, CVMNOZ 

, A8SAV 

, A9SA  V 

, A28SAV, A29SAV 

1 1 A0 

COMMON 

/ FRONT/ 

1 1 10 

IT  1 

,P  1 

.HI 

.SI 

,T2 

,P  2 

. H2  , S 2 

1 1 80 

3T2 1 

,P71 

, H?  1 

,S?1 

«T3 

. P3 

. H3  , S3 

1 1 90 

3T4 

»P4 

,H4 

,S4 

,T5 

,P5  ... 

, H5  . ,S5 

1200 

4T55 

,P  55 

,H55 

. S55 

, BLF 

.BtC 

, RI.PU  , BI.OR  . 

1 ?l  0 

5CNF 

, PRF 

, ETAF 

, WAFC 

,waf 

.NA3 

,NG4  ,FAP4  , 

1220 

6CNC 

,PRC 

, ETAC 

, WACC 

, MAC 

. ETAb 

, DPCOM  ,DUMF  , 

17  30 

7CNHP 

, FTATHP 

.DHTCHP 

,DHTC 

. BLHP 

. WG5 

,FA R5  ,CS 

1240 

8CN1.P 

.ETATLP 

.DHTCLP 

. DHTF 

, BFJ.P 

, WG55 

, FAP55  .HPFXT  , 

r>50 

OAM 

, * LTP 

, FT  AR 

,ZF 

,PCNF 

.ZC 

, PCNC  ,WFB 

1 260 

ATFFHP 

.TFFLP 

, PCBLF 

, PCBLC 

, PCBLDU 

.PCBLoB 

, PCBLHP , PCB LLP 

1 o70 

COMMON 

/ SIDE/ 

1 ?R0 

XXPI 

,XWAF 

, XWAC 

, XBLF 

, XB1.DU 

,XH3 

, DUMS1  .DIJMS2  , 

1 290 

XXT21 

,YP2  1 

, XH2  J 

. XS2  1 

,T23 

,P23 

,H23  ,S?3 

1 300 

3T24 

« P 24 

,H24 

.S24 

,T?5 

,P25 

. H25  .S25 

1310 

4T28 

.P28 

,H28 

,S2« 

, T29 

, P29 

, H29  ,S29 

1 370 

5WAD 

,WFD 

,WG24 

.FAR74" 

, FTAD 

.DPDUC 

. BYPASS, DUMS3  . 

1 330 

6TS28 

.PS28 

,V28 

.AM  28 

, TS29 

,PS29 

,V29  , AM29 

1 340 

COMMON 

/ RACK/ 

1 350 

XXT55 

, XP55 

, XH55 

. XS55 

.XT  25 

.XP2~ 

, XH  25  ,XS25  , 

1360 

XVWFB 

. XWG55 

, XFAR55 

.XWFD 

, XWG24 

,XFAR24 

.XXPl  .DUMB  , 

13  70 

3T6 

,P6 

,H6 

, S6 

,T7 

.PI 

,H7  ,S7 

1380 

4T8 

.P8  “ 

,H8 

.58 

,T9 

;P9'  " 

,H9  ’ ,S9  " . 

1 390 

5WG6 

,WFA 

,WG7 

.•FAR7 

. FTAA 

.DPAFT 

,V55  ,V25 

14  00 

6PS6 

,V6 

, AM6 

. TS7 

. PS7 

,V7 

, AM7  , AM25  , 

1410 

7TSR  ' 

,PS8 

,V8 

,AMR  " 

,TSO 

”,PS9 

;V9  ~ , AM9 

1420 

rva 

,FRD 

. VJD 

,FGMD 

, VJM 

, FGMM 

, FGPD  ,FGPM  , 

1430  __  9FGM  , FGP .JUFT  ,WGT  .FART  ,FG  _ tFN  ,SFC_ 

1440  COMMON  /DFRI  VX/2lU'X,Z2UX,  ZIUY,Z2ilY,  Z2UZ  I'.CPA TCPF'.CPFX.AKEX. 

1450  i REX,  HUFAR,  RUFA  R.  PUFAR,  ZFU,  PCNFU,  XT4 , ZCU , PANCU.TFFHPU , 

1 4A0_  ATFFF.PU 

1 470  DIMENSION- VAR( 6T,DEL(6),E RRBY6) , DEI  VAR<  6 ) ,EMAT (6,6),  WAT(6) 

1480  |AMAT(6),VU(6) 

1490_  _ J3ATA  A WORD/6HENGBAL/ 

1500  DATA  VDELTA.VLIM.VCHNGE.NOMIsX/ 

1510  1 0.001,0.100,0.850,4/ 

1520  _ CALL  INPUT 

1530  if( i Nit. eo. n gcTto  50 

1540  TFFHP-TFHPDS 

1 550  TFFI.P-TFLPDS 

"iS40#50  '*  L(X')PFR*0’  ‘ “ 
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Table  5.  Sample  Listing  of  Subroutine  ENGBAL  (with  Analytical 

Derivative  Logic)  (Subroutine  Not  Checked  Out)  (Continued). 

Ih7?l  NUMMAP-0 

I5R0  NOM I SS-0 

I 59  0#  I LOOP-0 

1600  MISMAT-0 

16  I td  NOMAP-0 

1620  I GO*? 

1630  DO  2 1*1,6 

1640  VMAT ( I )*0. 

1650  AMAT(I)*0. 

1660  VU(I)-0.0 

1670.,.  DP!.VAR(I)*0, 

1680  DO  ?L*I,6 

1690*2  EMATC I ,L) “0. 

1700*3  , LOO  PFR*  LOOP  ER-M  . . 

1710  CALL  COFAN 

1 720  WORD-AWORT 

1730  IF(L(X)PER.GT.  ITRYS)  Go  TO  20 

' 1 740  IFCNOMAP.CT.0)  GO  TO  I 

1750  NUMMAP-0 

1760  55  IDERV*0 

177 0 VU( 1 >=100.0*ZFU 

1780  IF( MODE. NF .3) VU( 2)*PCNFU 

1790  IF(MODE.EO.3)VU(?)®XT4*T4/ipi,0 

1800  VU(3)= I 00.0*ZCU 

1810  IF(MODE. NE. 1 )VU(4)-PCNCU 

1820  IFCMODE. EO. I )VU(4)«XT4*T4/10.0 

1830  ' VU(5)*TFFHPU 

1840  VIJ(6)=TFFLPU 

1850  _ IF(  IDERV.NE.  DGO  TO.  200  _ 

" i860  GAO.  >ERR 

1870  DO  200  1-1,6 

1880  ERR(  I )«ERRB(  I ) ♦ERRU(  I ) 

”1890'  GO  TO  210  

1900  200  IDERV-1 

1910  910  no  4 1-1,6 

1920  ' " IF(ABS<ERR(  I > .'GT.TOLALL ) GO  TO' 5 ' 

1990  4 CONTINUE 

1940  CALL  PERF  _ __  _ 

1950"  ' CALL'ERROR 

1960  5 IF(LOOP.GT.0)  GO  TO  .7 

1970  DO  6 1-1,6  

1980  ERRB(  I )*ERR(  I) 

1990  DEL (I)- 1.0 

2000 GO  TO  9 _ 

2010  7 DO  8 1-1,6 

2020  8 EMAT ( I ,LfX)P)— ( ERRB  ( I )-ERR(  I )) /DELI LOOP) 

2030  _ yU(L(K)P)-0.0 

2040"  9"  L(X)P— L (X)P+ 1 

2050  IF! L(X)P . GT .6 ) GO  TO  41 

2060 VU(L(X1P)«L.0 

' 20}0  10  2FU-VU(  I)/ 100.0 

2080  IF(MODF.NF.3)PCNFII— VU(  2) 

2090  _ IF(MODE.fO.3)XT4-VUi^)*i0s0/T4 

2100  ” ZCU-VU(3)/ 100,0 

2110  IF(MODE.NE. I JPCNCU-VU ( 4J 

2 120 IF! MODE.  EO.  1 )XT4-VUI4.)  *J0._0/T4 

2130'  " TFFHPU-VU(  57 
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Table  5.  Sample  Listing  of  Subroutine  ENGBAL  (with  Analytical 

Derivative  Logic)  (Subroutine  Not  Checked  Out)  (Continued). 

2140  TFFI.  PU*VIJ(  6) 

2150  IF(  ZF.I.T.  0. ) ZF-0.05 

2160  IF<ZC.LT.P.)  ZC-0.05 

?l 70  GO  TO  < 1 ,3) ,100 

2 I 80#  II  DO  12  I »|,6 

?l  90#  I 2 AMAKD  — FRRP(I) 

2? 00  DO  14  1*1,6 

2?  10  IZFRO-0 

2220  DO  13  LOOP*  1. 6 

22  30#  13  IF  ( EM  AT  (I  .LOOP)  . FQ.'0. ) I ZFRO-  I ZFRO+  1 

2240  IF( IZERO.LT. 6)  GO  TO  14 

2250  WRITEC6.  100)1 

2260  1.OOPFR-ITRYS+I00 

2270  GO  TO  20 

2280# 14  CONTINUF 

2290  DO  16  LOOP-1, 6 

2200  IZERO-0 

2310  DO  15  I-f,6  ■■■'  * ■ 

23  20#  15  IF(  EMAT(I  ,l.CK)P)  .EQ.0. ) I ZERO-I  ZERO*  1 

2230  IF.' IZFRO.LT.  6)  GO  TO  16 

2340  WRITE(  6,  .10 I )L(X>P  

22  80  !.<X)PfR»ITRYS-M00 

2360  GO  TO  20 

2370#  16  * CONTINUE"  - --  - 

2380#  17  CAM.  M ATR I X ( FMAT, VM AT , AMAT) 

2390  LBIG-0 

2400  VaRBTG*0.0  *“  ‘ “ ~ — - “ * - 

24  10  DO  18  L*  1,6 

2420  ABSVAR*AB5(VMAT(D) 

2430  " IF(ABSVAR^LE.VLIM*VAR(L))  GO  TO  18 * “ 

2440  IF(ABSVAR.LE.VARBIG)  GO  TO  18 

2450  LBIG-L 

2460'  ~ VARBIG-ABSVAR  “ ~ 

2470# 18  CONTINUF 

2480  VRATIO-1.0 

2490  ” IF(LBIG.GT.0)  VRATlO*VLIM*VARf  LBIGJ/'VARBIG  ‘ * 

2500  FRR AVF—0.0 

2510  VMTAVE-0.0  __ 

2520  DFLAVF-0.0 

25  30  DO  19  L*  1,6 

2540  _DF.I.VAR(l.>«VRATIO*VMAf(L) 

2550  ERRAVE«ERRAVE+ABS<AMAT<L‘>)  '*  " 

2560  VMTAVE*VMTAVE+ABS( VMAT(L) ) 

25  70  DFLAV&;I)ELAVe*-AB5(DFLVAR(L)) 

2580#  19  “ V A R < L ) • VA  Rf  L ) ♦DEL  V A R ( L ) ~ 

2590  FRRAVE-FRRAVF/6. 

2600  VMTA VE-VMTA VE/6. 

2610  "DEi.AVEiDELAVE/6." 

2620  IF(MISMAT.GT.0)  GO  TO  3? 

2630  _ IF(  NOMI SS.  EQ.0)  MISMAT-i! 

2640  I F f MTSM A T . FO. 0T  “ I G0*1 

2650#20  WRITE(8,102)  LOOPER 

2660  DO  21  I* 1,6 

2670# 2 T XRITEffl,  1 03rAMAT(I ) , (EUATHTLIVL-T ,67YVMATU ) ,DELVAR< H , VAR(I ) 

2680  WRITE! 8, 104)  ERRaVE, VMTaVE.DELAVE 

2690# 22  IF( LOOPER.LT. I TRYS)  GO  To  10 

— 2700  - * TAIL  “ERROR 
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Table  5.  Sample  Listing  of  Subroutine  ENGBAL  (with  Analytical 

Derivative  Logic)  (Subroutine  Not  Checked  Out)  (Concluded). 


2710  RETURN 

2720030  VKTAVX-VMTAVE 

2730  DC)  31  1-1,6 

2740031  AMATCI ) — ERR( I ) 

2750  GO  TO  17 

2760032  WRITE! 8 , 105 ) AMAT,ERRAVE, DELVAR.DFLAVE, VMAT.VMTAVE, VAR 

2770  MISMAT-M ISMAT+  1 

2780  IF!VMTAVE.LT.VCHNGE*VMTAVX)  GO  TO  22 

2790  WRITER,  106) 

2800  IF(MISMAT.LT.NOMISX)  .‘IOMISS-  I 

28)0  MISMAT-0 

28  20  LOOP-0 

2830  I GO- 2 

2840  GO  TO  55 

28500100  FORMAT!4H0ROW,  I 2,  I6H  IS  ZERO  IN  EMAT) 

28600)01  FORMAT! 7H0COLUMN, 12, I6H  IS  ZERO  IN  EMAT) 

28700102  FORMAT!  8HB  ERRB , 28X  23HE  RROR  MATRIX  A.FTER  LOOP . 14. 29X4HVMAT, 
2880  16X6HDELVAR,  7X I4HVARI  ABI.ES  $ $S  S$) 

28000103  FORMAT! IH0,F8.4,8X6FI0.4, I0XF1 0.4.F  1 1. 4 ,4XF 1 1 .4 ,6HS S$ $$$ ) 
29000  1 04  FORMAT!  IH0, F8.4.32X I4HAVERAGE  VALUES, 3 IX,?F  ll.4,6H$S$SSS) 

29100105  FORMAT! I 2H0- AMAT, 7F 16. 6, 6HS SS SSS , 

2920  I/,  I 2H DELVAR,7F|6.6,6HSSSS$S, 

2930  2/,  12H VMaT,.7F  16,6, 6HSSSSSS, 

2940  3/,  ' 1 2H  VAR.6F  1 6.  6, 6H  SSSSSS) 

29500  1 06  FORMAT!  1H0.50X22HCHANGE  T!X>  SMALLSSSSSS) 

2960  FND_ 
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A. 2 MAIN  COMPRESSOR  SUBROUTINE  DERIVATIVES 

A listing  of  the  main  compressor  subroutine  (COCOMP)  is  given  in  Table 
6.  A flowpath  showing  the  inputs  and  outputs  required  on  the  derivative 
path  is  shown  in  Figure  4.  As  can  be  seen  from  the  figure,  the  inputs  to  the 
subroutine  (COCOMP)  are  obtained  from  four  different  sources.  These  sources 
are : 

1)  Calculations  from  the  upstream  components.  These  calculations 
determine  the  derivatives  of  the  thermodynamic  properties,  i.e., 
dP21 , dT21 , dH21 , and  dS21 . 

2)  Inputs  from  the  control  logic  (subroutine  ENGBAL) . This  logic  will 
increment  the  independent  deck  variables. 

3)  Derivatives  calculated  from  the  compressor  map.  The  derivatives 
of  all  the  dependent  map  variables  (PRC,  WACC,  ETAC)  with  respect 
to  the  independent  variables  (ZC,  CNC) . 

4)  Values  of  thermodynamic  properties  and  derivative  values  returned 
from  the  thermo-subroutines. 

All  of  these  input  values  could  be  returned  either  through  expanded 
argument  lists  or  through  labeled  common.  The  outputs  from  the  subroutine 
include  the  derivatives  of  the  outlet  thermodynamic  properties  as  well  as  the 
flow  derivatives  (including  parasitic  flow  derivatives). 

The  procedures  for  differentiating  the  subroutine  will  now  be  given.  The 
source  equation  to  be  differentiated  will  be  given  first;  the  differentiaton 
will  then  follow  together  with  a FORTRAN  statement  of  the  derivative  formula. 

A listing  of  the  subroutine  (COCOMP)  with  analytic  derivatives  is  given  in 
Table  7. 

Note  that  in  the  equivalent  FORTRAN  statements  an  X prefix  has  been  used 
to  indicate  a log  derivative  (i.e.,  XT  ■ dT/T),  a trailing  U to  indicate  a 
differential  (i.e.,  HU“dH) , and  an  imbeded  U to  indicate  a derivative  (i.e., 
PRCUZC  - 3PR/CC). 


THETA  - SQRT  (T21/518.668) 

Equation 

dTHETA 

1 dT21 

Derivative 

THETA 

" 2 T21 

THETAU 

- 0.5  * THETA  * XT21 

FORTRAN  Statement 
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i 000 

10  I 0 
|0 70 
1030 
(1*40 

1050 

1060 

1070 

1080 

10O0 

1 1 00 
11(0 
1 1 20 
1 130 
1140 


Table  6.  Listing  of  Subroutine  COCOMP  (Main  Coapressor). 

SUBROUTINE  COCOMP 
COMMON  / All/ 

IWORD  , IDES  .JOES  , KDFS  .MOD*  .INIT  , IDIIMP  . I AMTP 
2IGASMX , IDBURN, IAFTBN, IDCO  . IMCD  . IDSHOC. IMSH0C.N07FLT 
3ITRYS  ,!.OOPFR, NOMAP  .NUMMAP.MAPFTXJ.TOf.A  1.1.,  FRR(  6) 

COMMON  /DESIGN/ 

IPCNFGU.PCNC0U.T40U  .DUMDI  .DUMP?  .DELFG  , DEI  FN  .DELSFC 
2ZFDS  .PCNFDS, PRFOS  , ETAFDS.NAFDS  ,PRFCF  .FTAFCF.WAFCF 
37.CDS  .PCNCDS, PRCDS  .ETACDS, WACDS  ,PRCCF  .FTACCF.WACCF 
4T40S  , WFBDS  , DTCODS . FTA8DS , HA 3CDS , DPCODS « PTCOCF , FTA BCF 

5TFHPDS.CNHPDS, ETHPDS.TFHPCF.CNHPCF.ETFPCF, rHHPCF,T2DS 
6TFLPDS.CNJ  PDS.FTLPDS.TFLPCF.CNT  PCF.FTTPCF.PHLPCF ,T?I DS 
7T24DS  ,WFDDS  , DT  DU  DS . ETA  DDS , HA  23DS . DPDU  DS , DTDUCF , FTA  DCF 
8T7DS  , NFADS  , DTAFDS , ETAADS . HG6CDS , DPAFDS , DTAFCF , ET a ACF 


9A55 


,A25  ,A6 


,A7 


,A8 


,A9 


,A?R  .A29 


1150 

APS55 

, AM55 

.CVDNOZ 

.CVMNOZ, A8SAV 

.A9SAV 

. A28SAV , A29SAV 

1 1 60 

COMMON 

/ FRONT/ 

11  70 

ITI 

.PI 

.HI 

.SI 

»T2 

,P2 

,H2  , S2 

1 1 80 

?T?I 

,P?1 

,H?1 

,S?1 

.T3 

.P3 

. H3  .S3 

1190 

3T4 

,P4  . 

,H4 

.S4.  _ 

. T5  . 

»P5 

,H5  ,S5 

1200 

4T55 

,P55 

,H55 

. S55 

• BLF 

,BLC 

, BI.DU  .BLOB 

1?  1 0 

5CNF 

,PRF 

,FTAF 

,WAFC 

.WAF 

.WA3 

,WG4  ,FAR4 

1220 

6CNC 

,PRC 

, ETAC 

. WAOC 

, MAC 

.ETAB 

, DP COM  ,DUMF 

1?30 

7CNHP 

.FTATHP 

.DHTCHP 

,DHTC 

,BLHP 

,WG5 

,FAR5  ,CS 

1240 

8CNLP 

.ETATLP 

, DHTCLP 

,DHTF 

, BLLP 

,NG 55 

, FAR55  .HPEXT 

1350 

9AM 

,AI.TP 

,ETAR 

..ZF  . 

,PCNF 

,zc 

.PCNC  ,WFB 

1 260 

ATFFHP 

.TFFLP 

, PCBLF 

.PCBLC 

. PCBLDU 

.PCBLOE 

.PCBLHP.PCBLLP 

1280 
"1290  ' 
1300 
13 1 0_ 
1320 
1330 
1340 
1 150 
1360 
1370 
1380 
1390 
1400 
'1410#  f 
)430#1 

1430 

1 440 

1450 

1460 

1470 

1480 

I 490# 100 
'1500 
1510#? 

15  20 

1530 

1540 

.1550  __ 

1560 


COMMON  / COMP/CNX( I5),PRX(15.15).MACX(I5.15),ETAX(15.15). 

I NCN, NPT(  15) 

DIMENSION  WLHt 2) * 

DATA  A WORT)  .W1.H/6HCOCOMP  ,6H  (I.O)  ,6H  (HI)  ✓ 

NORD-AWORD 

" THETA«S0RT(T2 1/5 TB. 6687  “ 

CNC- PCNC/ (100. *THETA ) 

IF(ZC.LT.0.)  ZC-0. 

IF<ZC.GT.i.)  Zc«l. ’ " 

CNCS-CNC 

_ CALL.  SEA RCH( ZC,CNC,PRC,HACCtETAC , 

' iCNXl  15.NCN.PRX' ( J.f ) >ACX(  lV I j . ETAX ( 7,'lT7NPT( T) ,~l 5," ' 15 , lGO> 
IF ( MODE. £0. 1 ) GO  TO  I 
IF  < (CNC-CNCS ) .GT.0. 0005*CNC>  MAPEDO  I 
IF(  iC^ErT.  T."0RriG(TrEQ.2)— WRITE!  8,  l000)”CNCS.>»!flYIOe)) 

000  FORMAT! 19H0*  * * CNC  OFF  MAP.FI0.4, 2XA6, 1 IH*  * *$$*$$$) 
MAC- WACC*P21 /THETA 

IFTI  NF.  H “OWTT)“? 


PRCCF-(PRCDS-I.) /(PRC-1.) 
FTACCF-ETACDS/ETAC 

tiaccf-wacds/wae 


WRITE!  6,  I00)PROCF,ETACCF.  WACCF.T2  IDS 
FORMAT! I 8H0COMPRESSOR  DESI0N.6X8H  PRCCF-.E15.8.8H  ETACCFj 

W 'tfACCFv  f iSTBVBff-  iBTSl * 

PRC-PRCCF*<  PRC- 1 . )♦ I . 

— 5T  A C*E T A CCF>ETAC 

MAC*  HAOCF*M AC 

CALL  THC()MP(PRC.ETAC.T21.H?I  .521  .P2J.T3.H3.S3.P3) 

-LEI  FCPLCL.GT.  0*JJUSCJ8LG*BAC 


WA3-NAC-BLC 


•jE.15.8. 
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Table  6.  Listing  of  Subroutine  COCOMP  (Main  Compressor)  (Concluded) . 


1570  RFJ)U-PCR?_DU*R!.C 

15R0  Bl.OB-PCBLOB*BLC 

1597  B!JdP*PCBI»P*R?.C 

1 600  BIXP»PCBU.P*BLC 

1610  IF( MODE. NE. I)  GO  TO  3 

16?0  IF(  ABS(CNC-CNCS)  .I.F.0.00UCNCS)  GO  TO  4 

1630  WRITE<8,2000)CNCS.CNC 

I640#?000  FORMATf 10H0CNC  WAS-  .PF5.fi,  FIH  AND  NOW-  .FI5.fi. 
1650  1 24H  CHECK  PCNC  INPUTSSSSSS) 

1660  CAM.  ERROR 

I 670#3  PCNC- I 00.*THETA*CNC 

I690#4  CAF.L  COCOMB 

1 690  RETURN 

1700 END 
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Table  7.  Sample  Listing  of  Subroutine  COCOMP  (with  Analytic 
Derivatives)  (Subroutine  Not  Checked  Out). 

1 000  SURROUTINF  CfCOMP 

If! IP  COMMON  / ALL/ 

10™  IWORD  .IDES  , JDFS  ,KDPS  ,MODF  .INIT  . IDUMP  , IAMTP 

1030  2IGASMX,  IDBURN,  IAFTBN,  IDCD  . IMCD  . I DSHOC.  I MSHOC . N07.FLT 

1040  3 1 TRY  S .LOOPFR. NOMAP  , NUMMAP.MAPEDG .TOLA LL, FRP(  A) 

1050)  COMMON  /DFSIGN/ 

1060  IPCNFGU.PCNC0U.T40U  , DUMD I , DUMP2  .DF1.FG  .DF1FN  .DEL RFC 

10170  2ZFDS  .PCNFDS.PRFDS  .FTAFDS.WAFPS  .PRFCF  .FTAFCF.WAFCF 

1080  3ZCDS  .PCNCDS, PRCDS  . ETACDS, WACDR  .PROCF  , FTACCF.WACCF 

1060  4T4DS  .WFRDS  ,nTConS.cTABnS,WA3CnS.nPa)n5,nTCOCF.CTAPCF 

I 1 00  5TFHPDS,CNHPDS,FTHPDS,TFHPCF, CNHPCF.FTHPCF.DHHPCF.T2DR 

I no  6TFI.PDS,CNLPDS,FTLPDS,TFLPCF, CNI.PCF, FTLPCF,  DHLPCF.T21 DS 

M 20  7T?4nS  .WFDDS  , DTDUDS , ETA DDS , NA  23DS, DPDUDR , DTDUCF , FTA DCF 

1130  9T7DS  .WFADS  .DTAFDS.ETAADS, WG6CDS,DPAr DS.DTAFCF .FTAACF 


1140 

0 A55 

.A  25 

. A6 

,A7 

, A8 

,A9 

,A2R 

.A  29 

1 150 

APS55 

, AM55 

.CVDNOZ 

, CVMNOZ 

, A8SAV 

, A9SA  V 

.A2PSAV 

, A29SAV 

1 1 60 

COMMON 

/ FRONT/ 

11  70' 

IT1 

,P1 

.HI 

.SI 

,T2 

.P2 

,H2 

,S2 

1 t 80 

2T21 

.P2I 

,H21 

, S2 1 

«T3 

,P3 

. H3 

.S3 

1 1 90 

3T4 

,P4 

«H4  _ 

• S4 

jT5  . 

.P5  . 

,H5 

, S5 

1200 

4T55 

,P55 

,H55 

. S55 

. BLF 

,BLC 

,B1  DU 

.BLOB 

12  10 

5CNF 

,PRF 

.FTAF 

.WAFC 

.NAF 

,WA3 

, NG4 

.FAR4 

1220 

6CNC 

,PRC 

, FTAC 

, WACC 

,HAC 

,ETAB 

, DPCOM 

.DUt'F 

16  30 

7CNHP 

.PTATHP 

.DHTCHP 

,DHTC 

, BLHP 

,NG5 

,FAR5 

.CS 

1 240 

BCNLP 

.ETATLP 

.DHTCLP 

, DHTF 

. BLLP 

,WG55 

, FAR55 

.HPFXT 

1 6S0 

9AM 

, ALTP 

.FT  AR 

. ZF 

,PCNF 

.zc 

, PCNC 

,HFB 

1 260 

ATFFHP 

.TFFLP 

, PCBI.F 

, PCBI.C 

, PCBLDU 

.PCBLOB 

, PCBLHP 

.PCBLLP 

1270  COMMON  / COMP/CNX(  I5),PRX(  15, 15) .NACX<t5, (5) ,ETAX( )5, )5) , 

I PRO  I NCN,  NPT<  15) 

1290"  COMMON  /DFRIVX/Z  IUX,Z21JX,ZIUY,Z2UY.Z2UZ  l,CP/~CPF'.CPFX,AKFX.' 

I 300  &RCX , HUFA  R , RUFAR , PUF  AR . ZFU, PCNFU , XT 4 . ZCU . PCNC1I.TFFHPU . 

13  10  &TFFLPU,XP2  1.XT2  1 .H2UJ.S2IU 

l 320 DIMFNS I'OfTWLHO)  ' 

1330  DATA  AW0RD.WLH/6HC0C0MP.6H  (LO)  ,6H  (HI)  / 

1340  IFUDERV.EO.  OGO  TO  3000 

1350  HoRD*A WORD  " ” 

1360  THETA»SQRT(T2 1/5  18.668) 

1370  CNC«PCNC/< 100.*THFTA> 

1380  iF(ZC.lT.'0n 'ZC-0.' 

1300  IF(ZC.GT.1.)  ZC*f. 

U00  _ cncs<;nc  __  

1 4 | 0 " CALL  SEA RCH( 2C ,CNC, PRC, HACC.ETAC , 

14  20  ICNX  ( I ) .NCN.PRXO  , I ) ,WACX<  I . U .ETAXX  I . I)  , NPT(  1 ) . 15,  15  . ICO) 

. 1 43  0 P RC.U  Zfi-Z  J UX 

1440  PRCUCN-Z1UY 

M50  WACCUZ-Z2UX 

1 460_  _ HACCU_N«Z2UY 

1470  ~ " FTACUZ*Z3UX 

14  80  FTACUN-Z3UY 

_U90 IF(^OpE-.EO,1)  GO  TO  J 

1500  I F ( ( CNC-CNCS ) . GT.0 . 0005*CNC)  MAPEDG-I 

1510  IF( IGO.E0.  I.OR. IG0.H0.2)  HRITE(8. 1000)  CNCS,HLH< 100) 

1520  _ FORMAT! 1OH0*  * * CNC  OFF  MAP.F1 0.4. 2XA6, 11H*  * *SSSSSS) 

1530  WAC»WAC&*P  2 I /THETA 

1540  IFUDES.'NF.  1)  GO  TO  2 

1550  PRCCF-CPRCDS-I .)  /(PRC-1 .) 

. -|560  FTACCF-FTACDS/ETAC 

38 


Table  7.  Sample  Listing  of  Subroutine  COCOMP  (with  Analytic 

Derivatives)  (Subroutine  Not  Checked  Out)  (Concluded). 

1570  WACCF“ WACD S/WAC 

1580  WRITC(6,  I00)PRCCF,ETACCF,WACCF,T2IDS 

1590  FORMAT(18H0COMPRESSOR  DESIGN, 6X8H  PRCCF«,  F 1 5. 8. 8H  FTACCF*, 

1600  1 P 1 5 . 8 , BH  W4CCF-.P15.8.BH  T21DS-, FI5.8 > 

1610  PRC-PRCCF*(PRC- I . )+1 

1620  ETAC-FTACCF*FTAC 

1630  WAC»WACCF*WAC 

1640  CALL  THC()MPCPRC.ETAC.T21.H21.S21.P2I.T3,H3.S3.P3> 

1650  CPA3-CPA 

1660  AKEX3-AKEX 

1670  REX 3 «R EX 

1680  IF( PCBLC.GT.0. ) BLC«PCBLC*WAC 

1690  WA3-WAC-BI.C 

1700  BLDU-PCBLDU*BLC 

1710  BLOB«PCBLOB*BLC 

1 7P0  BLHP-PCRLHP*RLC 

1730  BLLP-PCBU  P*BLC 

I74P)  IF(M()PEVNF.  1)  GO  TO  3 

1750  IF( ABS (CNC— CNCS) .LE. 0.001 *CNCS)  GO  TO  4 

1760  WRITE(B, ?000)CNCS,CNC 

1770  FORMAT ( 10H0CNC  WAS*  ,E15.8,.I1H  AND  NOW-  ,F15.R, 

1780  124H  CHECK  PCNC  INPUTSSSSSS) 

1790  CAM.  ERROR 

1 B00  PCNC- 1 00. *THETA *CNC  

18 10  GO  TO  4 

I B?0C  DERIVATIVE  PATH 

1830  3000  CONTINUE 

1540  THETAU=0.5*THETA*XT21 

1850  CNCU-CNSC*(PCNCU/PCNC— 0„ 5*XT21 ) 

I860  ' PRCU»PRCCF*( PRCUZC*ZCU+PRCUCN*CNCU)  ” 

IB70  ETACU-ETACCF*.(  NTACUZ*ZCU«-ETACUN*CNCU> 

1580  _ WACCU»WAOCF*(WACCUZ*ZCU*WACCUN*CNCU>  _ 

1 590  WACU-WACCU/WACC+XP2I -THETAU/THETA 

1900  BLCU-0.0 

1910 IE(PCBLC.GT.0)BLCU-PCBLC*WACU.. _ _ 

1920C  FLOW  DERIVATIVE 

1930  WA3U*WACU+BLCU 

1940C  _ PARASJTIC  FLOW  DERIVATIVES 

1950'  BLDUU-PCBLDU*BLCU 

I960  BL0BU-PCBJ.0B*BLCU 

1970  BLHPU-PCBI.HP  *BLCU 

1980  Bli  P U-PCB  LLP  *BLCU 

1990C  DERIVATIVES  OF  EXIT  THERMODYNAMIC  PROPERTIES 

2000  XP3-XP2J+PRCU/PRC 

20 1 0 XT 3-  XT 21  ♦ < AK FX3-  1 . 0_>  /AKE X3*PRCU/  ( PRC-EX  P ( 1 .0/AKEX3*  “ 

20  2 0 &ALOG (PRC ) ) )— ETACU/ETAC 

20  30  _H  31  I-  C P A 3*T  3*  XT  3._ 

2040  S3U-CPA3*Xf3-REX3*X~P3 

2050  4 CALL  COCOMB 
2060  RETURN 

2070 mi 
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2)  CNC  « PCNC/ ( 1100  * THETA) 


Equation 


dCNC  _ dPCNC  dTHETA 
CNC  “ PCNC  “ THETA 


dCNC  _ dPCNC  1 dT21 
CNC  “ PCNC  ' 2 T21 


Derivative 


CNCU  - CNC  * ( PCNCU/PCNC  - 0.5  * XT21) 


FORTRAN  Statement 


3)  The  following  compressor  map  derivatives  are  returned  from  SEARCH 
PRCUZC  - 3PRC/3ZC)cnc  - ZlUX 
PRCUCN  - 3PRC/3CNC)zc  - Z1UY 
WACCUZ  - 3WACC/3ZC)cnc  - Z2UX 
WACCUN  - 3WACC/3CNC)zc  - Z2UY 
ETACUZ  ■ 3ETAC/aZC)CNC  - Z3UX 
ETACUN  - 3ETAC/3CNC)Zc  - Z3UY 


4)  The  total  differentials  are  then  obtained  as  follows: 


dPRC  - 


3PRC\ 
aZC  /CNC 


dZC  + 


»CHCj 


zc 


dCNC 


PRCU  -PRCUZC  * ZCU  ♦ PRCUCN  * CNCU 


The  efficiency  and  flow  differentials  are  obtained  in  a similar  manner. 
ETACU  - ETACUZ  * ZCU  ♦ ETACUN  * CNCU 
WACCU  - WACCUZ  * ZCU  ♦ WACCUN  * CNCU 


5)  WAC  - WACC* P2 1 / THETA 

dWAC  _ dWACC  dP21  dTHETA  _ dWACC  dP21  Q 5 dT21 

WAC  ’ WACC  P21  " THETA  " WACC  P21  ~ * T21 

WACU  - WAC*(WACU/WACC  XP21  - 0.5*XT21) 
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6)  The  correction  (or  scale)  factors  must  now  be  applied  to  pressure  ratio, 
efficiency,  and  flow. 

PRCU  - PRCCF  * PRCU 

ETACU  - ETACCF  * ETACU 

WACU  - WACCF  * WACU 

7)  The  parasitic  flow  derivatives  are  obtained  as  follows: 
dBLC  - 0 

(BLCU  - 0) 

IF( PCBLC  .GT.  0)  BLCU  - PCBLC  * WACU 

WACU  - WACU  - BLCU 

BLDUU  - PCBLDU  * BLCU 

BLOBU  - PCBLOB  * BLCU 

BLHPU  - PCBLHP  * BLCU 

BLLPU  - PCBLLP  * BLCU 

8)  The  derivatives  of  the  compressor  exit  thermodynamic  properties  can  be 
obtained  in  the  following  manner. 

PRC  - P31/P21 

dP3  , d?21  dPRC 
P’-  P21  PRC 

XPj  » XP21  ♦ PRCU/ PRC 


dT3  _ dT21  + y d(PR)  _ dh^ 
T3  T21  PR-PRlA  He 
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ETACU 
ETAC 

dH3  = Cp  dT3 

H3U  = CPEX3  * T3  * XT 3 

dS3  = Cp  dT3  - RdP3 
T3  P3 

S3U  = CPEX3  * XT3  - REX3  * XP3 


(AKE3-1 )/ AKEX3  * PRCU 
PRC-EXP(1.0/AKEX3  * ALOG(PRC) 


4.3  MAIN  COMBUSTOR  SUBROUTINE 

A listing  of  the  main  combustor  subroutine  (COCOMP)  is  given  in  Table 
8,  A flowpath  showing  the  inputs  and  outputs  required  on  the  derivative 
path  its  shown  in  Figure  5.  As  was  the  case  with  the  main  compressor  inputs 
are  obtained  from  four  different  sources.  Note  that  the  combustor  inputs 
from  upstream  are  the  outputs  from  the  main  compressor  subroutine.  A listing 
of  the  subroutine  (COCOMP)  with  analytic  derivatives  is  given  in  Table  9, 

The  procedure  for  differentiating  the  subroutine  will  now  be  sunanarized. 
As  was  the  case  with  the  compressor  the  source  equation  will  be  given  first 
“cl lowed  by  the  derivative  formula  and  its  FORTRAN  equivalent: 

P3PSI  = 14.696  * P3  (Source  Equation) 

dP3PSI  = 14.696  * dP3  (Derivative) 

P3PSIU  = 14.696  * XP3  * P3  (Equivalent  FORTRAN  Equation) 

2.  WA3C  = WAC  * SQRT  (T3)/P3PSI 

dWA3C  _ dWAC  1 dT3  dP3PSI 
WA3CC  WAC  + 2 T3  “ P3PSI 

WA3CU  = WA3CC  * (WA3U/WA3  + 0.5*  XT3  - XP3) 


3.  DP COM  = DPCODS  * (WA3C/WA3CD) 
dDPCOM  - DPCODS/WA3CDS  * dWA3C 
DP CO MU  - DP00DS/WA3CDS  * WA3CU 


4.  IF(DPC0M  .GT.  1)  DP00M  - 1 

IF(J)PCOM  .EQ.  1)  DPOOMU  - 0.0 
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l c *00 
1010 
1020 
1030 

t 040 

1050 

1060 

1070 

I0R0 

1090 

1 I 00 

I MU 
1 1 90 
I I 30 

I 1 40 
1150 
I I 60 
1 1 70 
I 1 B0 
1 1 00 
I ?00 
I?  10 
I 220 
12  30 

1940 

1250 

19*0 

1 270 
1 980 

I 290 

1300 

1310 

1320 

1 330 

1340 

1350' 

1360 

1370 

1380 

I 390# 1 

1 400 

14  I 0 " 
1420 

1 4 30#  2 

|'440~  " 

| 4 50 
J_460#3 
■"'1470 
1480 
I490_ 
f500 
I510#4 
1890 

15  30#5 

1540 

1 550  _ 

1560 


Table  8.  Listing  of  Subroutine  COCOMB  (Main  Combustor). 

SUBROUTINE  COCOMR 
COMMON  / ALL/ 

IWORD  .IDES  ,JDES  , KDES  , MODF  .INIT  .IDUMP  , I AMTP  . 
2ICASMX, IDRURN, I AFTBN. IDCD  . IMCD  . I DSHOC. IMSHOC . NOZFLT. 
3ITRYS  .IOOPER, NOMAP  , NUMMAP.  MAPEDG.TOI.A  I.I.,  FPR(  6 ) 

COMMON  /DFSIGN/ 

1 PCNFGU .PCNCGU ,T4GU  .DUMDI  , DIIMD2  .DELFG  . DE1.FN  .DFLSFC. 
2ZFDS  .PCNFDS.PRFDS  .ETaFDS, WaFDS  .PRFCF  , FTAFCF . WAFCF  . 
3ZCDS  , PCNCDS, PRCDS  . ETACDS, WACDS  .PROOF  . FTACCF .WACCF  . 
4T4DS  , WFRDS  , DTCODS, ETABDS , WA3CDS .DPCODS, TTCOCF, ETARCF , 
5TFHPPS .CNHPnS, ^THPnS.TFHPCF.CNHPCF, FTHPCF.DHHPCF ,T2DS  . 
6TFLPDS.CNI  PDS,  FTLPDS.TFLPCF,  CNLPCF , ETI-PCF,  PHI  PCF.T2  IDS  , 
7T94DS  .WFDOS  .DTDUDS. FTAnDS, WAPTDS.DPDUDS.PTDUCF .RTADCF, 
RT7DS  , WFADS  , DTAFDS. ETAADS, MG6CDS, DPAFDS, DTAFCF, ETAACF , 


QA55 

, A 25 

,A6 

.*7 

. A8 

,A9 

, A28  . A 29 

APS55 

, AM55 

, CVDNOZ 

.CVMNOZ 

, A8SA V 

. A 9SA  V 

,A2PSAV,A29SAV 

COMMON 

iTl 

/ FRONT/ 

.Pi  .HI 

.51 

,T2 

.P? 

,H2  .S2 

2T2I 

,P  2 1 

,H2  1 

, S2  1 

,T3 

,P3 

,H3  .S3 

3T4 

,P4 

,H4 

.S4 

,T5 

• P5 

• H5  ,S5 

4T55 

,P55 

, H55 

, S55 

, BLF  ‘ 

,blc 

, Bl.DU  , BLOB 

5 CNF 

,PRF 

, ETAF 

. WAFC 

, WAF 

,W  A3 

, WG4  ,FAR4 

6CNC 

, PRC 

,ETAC 

. WACC 

, MAC 

.ETAB 

.dpcom  ,dumf 

7CNHP 

.FTATHP 

.DHTCHP 

,DHTC 

, BLHP 

,HG5 

, FAP5  ,CS 

RCNLP 

.FTATLP 

.DHTCLP 

,DHTF 

,BLLP 

,WG55 

, FAR55  .HPFXT 

9AM 

, ALTP 

, FTAR 

.ZF 

,PCNF 

,zc 

, PCNC  , WFB 

ATFFHP 

.TFFLP 

, PCRJ.F 

, PCRLC 

, pcrldu 

.PCBLOR 

« PCRLHP  , PCB 1 LP 

COMMON  / COMB/PSI ( 15) ,DELT (15, 15),ETA( 15, 15), NPS, NPT (15) 
DIMENSION  0(9) ,DUMBO( 15,15)_ 

DATA  AW0RD/6HC0C0MB/  ' ~ ' 

HORD-AWORD 


O<2>=0. 

O(3)«0. 

P3PS 1= 14 ,696*P3 
NA3C»WA3*S0RT(T3)/P3PSI 
IFf IDES.EQ. 1 ) NA3CDS-HA3C 
DPCOM»DPC()DS*  ( WA3C/WA3CDSJ 
IFfDPCOM.GT.  1 . ) DPCOM^J 

P4*P3*( 1 . -DP COM) 

I Ft T4.GT. 3999. ) T4«3999. 

IF(T4.GF.  1000. ) GO  TO  ? __ 

T4»1 000. 

IFf MODE. EO. 1)  MAPEDG* 1 
DTCO-T4-T3 

'ifudes.ne:  ir&O  TO  3“  ‘ 


DTCOCF-DTCODS/DTCO 
DTCO-DTCttCF *DTCO 

P 3P  S IN«P3P  SI  ' " "*  ~ 

CALL  SEARCH ( -1 . .P3PSIN.DTC0.ETAB. DUMMY. 

IPSIf  D.NPS.DELTf  I.D.ETAf  I,  1),DUMB0(  1,1  ),NPT(  0.15, 15. 1 GO) 
iFd'GO.  FO.  7)  CAtr  ERRfiR  "" 


IFf  IDES.NE.  I)  GO  It)  5 
FTABCF-FTABDS/CTAB 

~ mB-ETABCFAETA'B' 

HV«(  (((!>(-. 4594317F-19*T4)-. 2034 1 1 6E- 15 ) *T4+*. 2783643F- 1 1 ) *T4 
I*.  205  I50IE^07J#T4^.  24531 16E“03)#T4-.9433296E— 01  )-#T4*.  IB45537E+05 
CAlL  “THERM0fT>4,TfA  sT4,  jfjf  i;*XX>V0,0.'0,0T 
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Table  8. 

Listing  of  Subroutine  COCOMB  (Main  Combustor)  (Concluded). 

1570 

FAR4* (HA-H3)/(HV*ETAB> 

15  80 

I F(F  AR4. 1.T  .0.  ) FAR4-0. 

1 590 

WFBX»FAR4*WA3 

1 AO0 

IF(  MODF'.NF.  ?)  GO  TO  8 

1610 

F RRW  * ( WF  B-  fcF  BX ) /WFB 

1 A70 

niR*SQRT(WFB/WFBX) 

1630 

CAM  A FOUI  R(f)(  1 ) ,T4 , ERRW, 0 . , 20.  . 0. 0001 . DIR.T4T,  I GO ) 

1640 

GO  TO  (6,9,7) , I GO 

1650#6 

T4-T4T 

1 660 

GO  TO  1 

1 670#7 

CAM  ERROR 

I680#8 

WFB*WFBX 

I690#9 

CAM  THERM0(P4,H4,T4,S4,  XX?,  1 ,FAR4,0) 

1 700 

WC4-WFB-HKA3 

17  10 

IF( IDES.FO. 1)  WRITE(6, 100)  WA3CDS, FTABC F,DTCOCF 

1 7?0# 1 00 

FORMAT( 17H0COMBUSTOR  DFSIGN.7XBH  WA3CDS-, FI  5. 8. 8H  FT ABCF*. F 15.8. 

1 730 

I8H  DTC()CF*,E  15.8) 

1740 

CAM  COHPTB 

1750 

RETURN 

1760 

_END 

— 

— ...  - 

ihispm*  is  best  ,u*im  raAciia«i 
from  OOP*  WRUISHED  TO  DQQ 


Table  9.  Sample  Listing  of  Subroutine!  COCOHB  (with  Analytic  Derivatives) 
(Subroutine  Not  Checked  Out). 


1000  SUBROUTINF  COCOMB 

IOIO  COMMON  / A!.!./ 

l*’0  IWOR0  ,IDFS  , JOFS  .KDFS  .MOD*  .INIT  . IPUMP  , I AMTP  . 

1030  2IGASMX.IDBURN, IAFTBN, IDCD  , IMCD  , IPSHOC, IMSHOC.NOZFLT. 

10^0  3 1 TRY  S .T.OOPFR.NOMAP  . NUMM* P.MAPCDG.TOLA  1.1  , FRR(  4) 

1350  COMMON  /DESIGN/ 

1060  IPCNFGU  .PCNCGU.T4GU  , DUMB  I ,DUMI)2  .DEI.FG  .CEIFN  .DE1.SFC. 

1070  2ZFPS  , PCN FDS, PRFDS  . FTAFDS, NAFDS  , PRFCF  . FTAFCF . WAFCF  . 

I05O  3ZCDS  , PCNCDS, PRCDS  , ETACDS, WACDS  ,PROCF  , FTACCF, WACCF  . 

1000  aTADS  , WFRPS  .DTCODS, FTABPS, NA3CPS,nPO)DS,rTC0CF , FTA RCF. 

1100  5TFHPDS.CNHPDS, ETHPDS, TFHPCF, CNHPCF.ETHPCF, DHHPCF.T2DS  , 

IM0  4TFLP  DS  ,C  NI.PDS , FTLPDS, TFLPCF, CNLPCF . FTLPCF, PHI PCF  ,T2 1 PS  . 

M?0  7TP4DS  .WFDDS  , DTDUDS.  ETADDS,  WA23DS.DPDUDS,  DTDUCF,  FTAPCF , 

1130  8T7DS  , WFADS  .DTAFDS.ETAADS.HGaCDS.DPAFDS.DTAFCF.ETAACF, 


1 i 40 

9 A55 

,A25 

,A6 

. A7 

,AR 

, A9 

,A2R  .A  29 

1 1 50 

APS55 

, A M55 

.CVDNOZ 

, CVMNOZ 

, A8SAV 

,A9SA  V 

• A 28SAV , A29SAV 

1 1 60 

COMMON 

/ FRONT/ 

1 1 70 

IT! 

.PI 

.HI 

.SI 

,T2 

,P2 

,H2  , S2 

1 ! 80 

2T21 

.P  2 1 

,H21 

. S2 1 

,T3 

.P3 

,H3  ,S3 

1 190 

3T4 

• P4 

.H4  ... 

,S4  ._ 

.15  . 

.P5 

,H5  ..  . S5. 

1200 

4T55 

,P55 

,H55 

. S55 

, BLF 

,BLC 

, BI.PU  .BLOB 

1^10 

5CNF 

, PRF 

,ETAF 

• WAFC 

, WAF 

• HA3 

.NG4  ,FAR4 

1220 

6CNC 

,PRC 

, ETAC 

, WAOC 

, HAC 

, ETAB 

, DPCOM  , DUMF 

12  ?0 

7CNHP 

.FTATHP 

.DHTCHP 

.DHTC 

,BLHP 

,WG5 

.FAR5  ,CS 

1 240 

RCNLP 

.ETATLP 

, DHTCLP 

,DHTF 

, BLI.P 

,WG55 

, FAR55  .HPEXT 

12  50 

9AM 

, ALTP 

,FTAR 

,ZF.  . 

,PCNF 

,ZC 

, PCNC  ,WFB 

1 260 

ATFFHP 

.TFFLP 

. PCBLF 

. PCBLC 

, PCBLDU 

.PCBLOB 

. PCBLHP .PCBLLP 

12  70  COMMON  / COMB/PSH  15),DELT(|5, 15).FTA<15,15),NPS,NPT( 15) 

1 2P0  COMMON  /PFR!  VX/Z1UX,  Z21JX,  Z1UY,Z  2UY, Z?UZ 1 ,CPA ,CPF,CPFX . AKFX , 

1290  & RFX, HUFA R,  RUFA R,  PUFAR.ZFU,  PCNFU,  XT4 , ZCU ,PCNCU, TFFHPU . 

1300  &TFFLPU,XP21,XT2J,H2UI,S21U,XP3,XT3,H3U.S3U,WA3H.WFBlJ 

13  10  DIMENSION  Q(9),DUMB0(I5, 15) 

13  20  DATA'  AWORP/AHCOCOMB/ 

1330  IFCIDERV.EQ. 1)00  TO  )50 

1340  WORD-ANORD 

‘1350  Q(  2-)  *0.  '■  ' “ 

1360  O(3)»0. 

1370  P3PSI*I4.696*P3 

13R0  WA3C-WA3*S0RT(T3)/P3P'sr 

1390  IFdDES.EO.  1)  HA3CDS-NA3C 

1400  DPC0M*DPC0DS*(HA3C/WA3CDS) 

"14  10  "'  IFFDPCOM.'GT.  I.)  DPCOM-:i: 

1 4 ?0  P4*P3*( 1,-DPCOM) 

_1430  1 I F( T4 . GT. 3999_. ) T4-3999.  

1440  IF(T4.GF.  J000.1  GO  T6~2 

1450  T4-1000. 

1460  IF(MODE.J0.i)_MAPEDG-L 

14  70  2 DTCO*T4—T3 

1480  IFXIDES.NE.  I)  GO  It)  3 

_1 4 90 DXCoCF-DTCODS/BICQ 

1500  3 DTCO-DTCOCF  *OTCO 
1510  P3PSIN-P3PSI 


1520  _ CALL  SEARCH(_-1  .,P3PSIN,DTC0,ETAB, DUMMY. 

1530  i ps i x i) ,nps,deLT(  i,  d;eta(  i,  n ,'dOmboc  i,i  r.jjptd),  15,1 '57160) 

1540  FTABUP»Z2UY 

1550  ETABUD-Z2UZ 1 

1 540  ‘ IF(  rG0.F&.7)  CALL  FRROR 
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Table  9.  Sample  Listing  of  Subroutine  COCOMB  (with  Analytic  Derivatives) 
(Subroutine  Not  Checked  Out)  (Continued). 


1 S7C*  4 IF<  IDFS.NF.  1)GO  TO  5 

1 5 SO  ETABCF*ETABDS/ETAB 

1590  5 FTAB=ETABCF*ETAB 

I *00  HV4*  ( ( (( <<-.45943l7F-l9*T4)-.20341 16E-I5 )*T4+.27R3643F-  t 1 >*T4 

15  10  I*. 205  150  IE-07) *T4-. 24  5 3 116E-03)>T4-. 94  3 3 296E-0 1 )*T4*. I845537F+05 

1670  CAM.  THFRPO(P4.HA,T4.XX1.XX2,0.0.0.0) 

1630  HA4*HA 

1640  FAR4* (HA-H3)  /( HV-*ETAB) 

1650  1F(F *R4. 1.T.0. ) FAR 4*0. 

1660  WFBX*FAR4*WA3 

1670  IF( MODF.NF.?)  GO  TO  B 

16B0  ERRW=(WFB-WFBX)/WFB 

1690  niR-SORT(WFB/WFBX) 

I7O0  CALL  AFOUIR(D( l),T4,ERRM,0.,20..0.0O0l*DJR.T4T.IGO) 

1710  GO  TO  (6,9,7) , IGO 

1770  6 T4*T4T 

1730  GO  TO  I 

1740  7 CALL  ERROR 

1750  8 WFB*WFBX 

17609  CALL  THERM0(P4,H4,T4,S4,XX2, I.FAR4.0)  _ 

1770  WG4-HFB+WA3 

I7R0  HUFAR4-HUFAR 

1790  RUFA R4*RUFAR 

1500  PUFAR4.PUFAR 

1510  REX 4* REX 

1820  CPEX4-CPEX  _ 

1830  CPA4*CPA 

1840  I F ( IDES.EO. 1 ) WRITE(6, 100)  WA3CDS.ETABCF .DTCOCF 

1950  100  FORMAT(  1.7H0COMBUSTOR  DESIGN, 7X8H  NA3CD5*,  E15'.8,8H  PTABCF*,  _ 

1860  ' IEI5.R.8H  DTCOCF-,E 15^81 

1870  GO  TO  200 

1RR0C  DERIVATIVE  PATH  _ _ _ 

1 890  150  P3PSiU*l4. 696*XP3*P3 

1900  WA  3C  U*  NA  3C*( WA  3U/H  A3+  0 . 5*  XT3-XP3) 

1910  DPC0MU*DPC0M*DPC0M*HA3CU/»»A3C  _ 

‘ 1920'  IFIDPCOM.EO. 1 )DPC0*U«0.0 

1930C  EXIT  PRESSURE  DERIVATIVE 

1940  _ XP4*XP3-DPC0MU/1  1S0-DPC(>M> 

1950  IF(T4.Eb.3999.)XT4-0.0 

I960  IF(T4.EQ. 1 000. ) XT4*0. 0 

1970  DTC0U*DTC0CF*(T4*XT4-I3#XTi3J[ _ __ 

1 9 ROC  EFFICIENCY  DERIVATIVE 

1 990  ETABU-ETABCF  *f  ETABUPAP3PSI U+ETABUD^DTCOU) 

20OOC EUEL..H.EATI  N.  VALUE. DERIVATIVE _ . 

2010  HVUT4-U  ( ( (-0. 2756590E—  18)  *T4-0.  J0 I7058E— 14)4T4 

20  20  4+0.  1 1 13457E-n)*T4+0;6.l54503E-07)*T4-O.4906?32E-03>*T4 

2030 Ar0._9433296E-01  

2040C  FUEL  TO  AIR  RATIO  DERIVATIVE 

2050  FAR4U-FA  R4* ( (CPA4/I HA4-H37-HVUT4/HV4)*T4*XT4- 

2060 AH3U/(HA4-H3)-ETAB)U/ETAB 

7070  if(faR4:eo-.0)far'4U-0;0 

2080  IF(M0DE.NE.2)G0  TO  170 

_ p090 FAR4U*— FAR4*WA3U/WA3  _ _ 

2100  XT4*T4*iH3U/(HA4— H37+ETABU/ETAB)  /(CPA4/<HA4»H33— HVUT4/HV4J 

2.110  170  CONTINUE 

2120  HFBU-FAR4*NA3U+*(A3*FAR4U 

2130'  HG4U«hEB  U4HA3U  “ " 
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Table  9.  Sample  Listing  of  Subroutine  COCOMB  (with  Analytic  Derivatives) 
(Subroutine  Not  Checked  Out)  (Concluded). 


2140C  DERIVATIVE  C)F  THE  EXIT  THERMODYNAMIC  PROPERTIES 

? , sei  H4U»CPH(4*T4*XT4+HUFAR4*FAR41J 

2 160  S4U«CPEX*XT4/<  1 .0-*-FAR4)-REX44XP4f(puFAR4-RIlFAR44 

2170  &AI.OG  (P4)  )*FAR4U 
21R0  200  CALL  COHPTB 
2190  RETURN 

FND 


5.  P4  - P3  * (1  - DPCOM) 

dP4  _ dP3  - dDPCOM 
P4  " P3  1 - DPCOM 

XP4  - XP3  - DP  CO  MU/  ( 1 - DPCOM) 


6.  IF(T4  .GT.  3999.)  T4  - 3999. 

IF(T4  .EQ.  3999.)  XT4  - 0.0 

7.  IF(T4  .GE.  1000.)  GO  TO  2 
T4  - 1000. 

IF(T4  .EQ.  1000.)  XT4  - 0 

8.  DTCO  - (T4  - T3)  * DTCOCF 
dDTCO  - (dT4  - dT3)  * DTCOCF 

DTCOU  « DTCOCF  * (T4  * XT4  - T3  * XT3) 


9.  The  following  combustor  map  derivatives  are  returned  from  subroutine 
SEARCH 

ETABUP  - 3ETAB/8P3PSIN)DTCO 
ETABUD  - 3ETAB/3DTCO)p3pslN 


dETAB  - ) dP3PSIN  ♦ ) 

*P3PSIN  DTC0  #DTC0  P3PSIN 


F.TABU  - ETABUP  * P3PSIU  + ETABUD  * DTCOU 


dDTCO 


10.  ETAB  - ETABCF  * ETAB 
ETABU  - ETABCF  * ETABU 


11.  HV4  - -0.459431 7E-1 9 * T46  - 0. 2034116E-15  * T45 
♦0.2783643E-11  * T44  ♦ 0. 2051501E-07  * T43 
-0. 24531 16E-03  * T42  - 0.9433296E-01  * T4 
♦0. 1845537E>05 
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HVUTA  = ( ( ( ( ( -0 . 2756590E-18)  * TA  - 0. 101 /058E-1A)  * TA 
+0. 1 11 3A57E-1 1 ) * TA  + 0 . 61 5A503E-07)  * TA 
-0.A906232E-03)  * TA  - -0 . 9A33296E-01 


12.  FARA  = (HA  - H3)/(HVA  * ETAB) 

dFARA  _ d(HAA-H3)  d(HVA  * ETAB) 

FARA  HAA-H3  ' HVA  * ETAB 

dFARA  _ dHAA  dH3  HVA  * dETAB  ETAB  * dHVA 

FARA  * HAA-H3  HAA-H3  " HVA  * ETAB  “ HVA  * ETAB 

dFARA  _ CPAA  dTA  dH3  dETAB  dHVA  dTA 

FARA  ' (HAA-H3)  “ (HAA-H3)  ' ETAB  _ dTA  HVA 

FARAU  = FARA  * ( CPAA/ (HAA-H3)  - HVUTA/HVA)  * TA  * XTA 

-H3U/(HAA  - H3)  - ETABU/ETAB 

IF( FARA  .EQ.  0)  FARAU  = 0 

13.  IF (MODE  .NE.  2)  GO  TO  170 
FARA  = WFB/WA3  (Constant  WFB) 

dFARA  = _ dWA3  r CPAA  dTA  dH3  dETAB  dHVA  dTA 1 

FARA  ” WA3  [_  (HAA-H3)  " (HAA-H3)  " ETAB  “ dTA  HvAJ 

FARAU  = - FARA  * WA3U/WA3 

dTA  = f d,H3  + d£TAB  - dWA3^  / ( CPAA  _ dHVA.  1 ^ 

^(HAA-H3)  ETAB  WA3  )/  \^HAA-H3)  dTA  HVA/ 

; 

XTA  - TA*(H3U/ (HAA-H3)+ETABU/ETAB-WA3U/WA3) /(CPA4/ (HAA-H3) -HVUTA/HVA) 

1A.  170  CONTINUE 

WFB  = WFBX  - FARA  * WA3 

dWFB  = FARA  * dWA3  + WA3  * dFARA 

WFBU  « FARA  * WA3U  + WA3  * FARAU 
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15.  WG4  - WFB  + WA3 

dWG4  - dWFB  + dWA3 
WG4U  - WFBU  ♦ WA3U 


16.  The  derivatives  of  the  tnermodynamic  properties  of  the  exit  are  obtained 
as  follows: 

a)  Enthalpy  Derivative 
H - (Ha  ♦ fHf )/(l  ♦ f) 


3H  \ Hf  - H 
'P,T 


8f/o  1 ♦ f 


HUFAR 


dH 

» V, 

dT 

*,  f 

♦ 2* 
3P 

X,f 

. _ 3H> 

df>  + 3f> 

1 

P,T 

dH 

- _L 

dHa 

f 

dHf 

+eHf-H 

df 

dT 

i+f 

dT 

l+f 

dT 

l+f 

‘ dT" 

dH 

CPA 

+ f 

rpf 

Hf- 

H df 

dT 

" l+f 

l+f 

Lrl 

+ T+f 

dT 

dH 

. CPA 

+ f*CPf  A 

Hf-H 

df 

dT 

l+f 

l+f 

dT 

dH  - CpdT 


<&) 


df 


H4U  - CPEX4  * T4  * XT4  + HUFAR4  * FAR4U 

where  HUFAR  and  CPEX4  are  returned  from  the  thermo  subroutine  PROCOM. 


b)  Gas  Constant  Derivative 
AMW  - 28.97  - 0.946186  * FARX 


REX  - i . 986375/AMW  - 1.986375/(28.97-0.946186  * FAR4) 


The  above  equations  define  the  molecular  weight  and  gas  constant  and 
appear  in  subroutine  PROCOM. 
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0.946186*R£X 

AMW 


3R£X\ 

3f  K,  T 


RUFAR4  - 

8f  ' P,T 

where  RUFAR4  is  returned  from  thermo  subroutine  PROCOH 


c)  Entropy  Derivative 
s - ♦ - R * ALOG(P)  - S(P,T, f ) 

where  « - »♦  C T , f ) 


l + f 


PUFAR 


& 

d*  - It)  dT  * It)  d£ 

f T 

d*  , d»a>fd»f  + 3*\ 

9 l+f  it) 

T 

But  d+  ■ cpdT/T  (by  definition) 


d* 


cPa  1 fcPf  ♦ i±\  df 


l + f 


8f/- 


d*  -C^-  dI-  +3^  df 
* l+f  T if  J 

'T 


dS  - d+  - - ALOG(P)  df 

p if 

is ' & f * It)  d£  - *¥ ' li AU>G(P)  d£ 

T 

S4U  - CPEX4  * XT4/U+FAR4)  ♦ PUFAR4  * FAR4U 

- REX4  * XP4  - RUFAR4  * ALOG(P4)  * FAR4U 
where  PUFAR4  is  calculated  in  the  thermo  subroutine  PROCOM 
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5.0  COST  ADVANTAGES  OF  THE  ANALYTICAL  DERIVATIVE  METHOD 


5.1  SELECTION  OF  FLIGHT  ENVELOPE  AND  MATRIX  OF  TEST  POINTS 


The  engine  data  matrix  used  in  the  current  Analytical  Derivative  pro- 
gram is  shown  in  Figure  6.  The  data  matrix  is  a modification  of  a similar 
matrix  found  to  be  useful  in  a typical  military  customer  deck.  It  contains 
a total  of  411  operating  points.  The  figure  shows  the  flight  envelope  with 
the  matrix  of  data  superimposed.  Each  of  the  locations  indicted  by  a circle 
consists  of  seven  points  while  those  identified  by  a triangle  consist 
of  five  points.  All  of  the  designated  locations  were  run  at  the  follow- 
ing power  levels: 

(1)  Maximum  Afterburner  Power 

(2)  Roughly  50  Percent  Afterburner  Fuel  Flow 

(3)  Intermediate  (Max  Dry)  Power 

(4)  Idle 

In  addition,  the  locations  indicated  by  triangles  were  run  at  one  part  power 
level  roughly  equally-spaced  between  Idle  and  Intermediate  (maximum  dry) 
power.  At  the  locations  designated  by  circles,  three  Intermediate  part 
power  levels  were  run  so  as  to  divide  the  range  between  Idle  and  Intermed- 
iate power  into  four  approximately  equally  spaced  intervals. 

All  points  were  run  at  standard  day  conditions,  and  the  ram  recovery 
was  in  accordance  with  MIL  Spec  5008B. 


5.2  COST  COMPARISON  FOR  INTERNAL  DECKS 


Computation  cost  comparisons  were  made  between  runs  using  analytical 
derivatives  and  runs  using  a reference  deck  which  had  not  been  modified  for 
analytical  derivatives.  Valid  comparisons  could  not  be  made  using  the 
same  deck,  even  though  numerical  derivative  capability  was  retained  in  the 
modified  deck,  since  the  computation  and  storage  of  partial  derivatives  in 
the  low-level  subroutines  extended  the  time  for  computing  base  passes  signifi- 
cantly. The  comparisons  were  made  using  the  411-point  matrix  of  flight  con- 
ditions, with  both  the  primary  (Engine  1)  and  alternate  (Engine  2)  engine 
models. 

Processor  time  was  measured  over  an  interval  that  included  all  the 
actual  engine  cycle  calculations,  but  not  the  program  initialization  or 
input-output  processing.  The  analytical  derivative  model  required  62  per- 
cent less  time  than  the  numerical  derivative  model  using  Engine  1,  and 
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61  percent  less  time  using  the  Engine  2.  Processor  time  is  not  a complete 
measure  of  the  relative  total  cost  of  running  the  models,  as  the  analytical 
derivative  model  requires  25%  more  computer  memory. 

The  majority  of  the  time  gain  (54%)  was  directly  attributable  to  the 
inclusion  of  the  analytic  derivatives.  About  8%  was  due  to  the  rewriting 
and  inclusion  of  a new  and  faster  version  of  the  subroutine  COFFEX.  This 
was  a relatively  new  subroutine,  and  the  coding  had  not  been  optimized. 

The  remaining  subroutines  were  not  new  and  the  coding  had  already  been 
opt imized . 

A total  cost  comparison  was  made  between  the  two  models,  each  run  from 
loader  output  files  without  any  subroutine  recompilation,  and  using  Engine  2. 
The  total  cost  of  the  run  using  the  analytical  derivative  deck  was  44% 
lower  than  the  run  using  the  numerical  derivative  deck  (38%  due  to  ana- 
lytic derivatives  and  6%  due  to  the  new  subroutine  COFFEX.  The  total  cost 
includes  not  only  the  processor  time  and  memory  residence  charges  associated 
with  performing  the  actual  cycle  calculations,  but  also  the  model  initiali- 
zation, input  processing,  and  output  formatting  and  conversion  charges.  The 
output  processing  cost  alone  has  been  measured  to  be  30%  of  the  total 
cost  of  running  the  analytical  derivative  deck. 

It  is  estimated  that  the  first  method  (i.e.,  the  complete  set  of  ana- 
lytic derivatives)  if  completely  implemented  internally  would  result  in  an 
(AEG  Evendale)  computations  cost  savings  of  about  $100,000  annually.  How- 
ever, the  many  second  order  effects  which  are  included  in  production  cycle 
decks  result  in  nonstandard  logic  in  the  higher  level  subroutines.  For 
this  reason,  about  $10,000  per  deck  is  estimated  for  inclusion  of  the  com- 
plete derivative  set.  In  addition,  continued  modification  of  component 
logic  would  require  continual  updating  of  derivatives  and  a cost  gain  would 
be  hard  to  realize. 

The  second  method  (low  level  derivative)  would  result  in  a cost  savings 
of  approximately  $75, 000/year . Moreover,  it  would  require  only  about  $500 
in  total  cost  per  deck  to  implement  the  method.  Method  2 Is  the  one  of 
choice  for  all  but  the  most  heavily  used  parametric  customer  decks. 


5.3  COST  COMPARISON  FOR  EXTERNAL  (CUSTOMER)  DECKS 


A total  cost  comparison  between  the  base  deck  and  the  deck  with  full 
analytical  derivatives  included  was  made  at  Wright  Patterson  Air  Force  Base 
on  a CDC  6600/CYBER- 74  computer.  Both  Engine  1 and  Engine  2 were  compared. 
Engine  1 showed  a 51.6%  savings  in  cost  while  Engine  2 showed  a 52.2% 
savings.  If  these  numbers  are  taken  as  fairly  typical,  a cost  savings 
of  about  52%  should  result  from  the  inclusion  of  the  method  in  an 
external  customer  deck.  The  difference  in  cost  savings  between  the 
external  and  internal  decks  can  be  attributed  to  the  greater  number  of 
output  parameters  printed  out  with  Internal  decks. 
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As  described  in  the  previous  section,  the  majority  of  this  gain 
(i.e.,  approximately  45%)  is  due  to  analytic  derivatives.  The  remaining 
6%  is  due  to  the  inclusion  of  the  new  faster  version  of  the  subroutine 
COFFEX. 


i 
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6.0  ASSESSMENT  OF  THE  GENERALITY  AND  UTILITY  OF  THE  APPROACH 


The  analytic  derivative  method  can  be  applied  to  any  cycle  deck  which 
uses  numerical  derivatives.  Based  upon  the  experience  at  AEG,  the  best 
approach  is  to  differentiate  that  portion  of  the  deck  for  vhch  the  coding 
remains  constant;  for  example,  the  subroutines  used  to  evaluate  thermo- 
dynamic properties  are  usually  good  candidates  since  their  structure  is 
fairly  constant  and  they  are  called  frequently.  This  approach  should  re- 
sult in  the  greatest  savings  per  unit  cost. 

For  the  Air  Force  Deck  (SMOTE)  the  deck  coding  appears  to  be  relatively 
constant  (based  upon  a comparison  of  the  current  coding  and  that  given  in 
Reference  3).  This  is  due  largely  to  the  neglect  of  second  order  effects 
such  as  tip  clearances,  etc.  Since  this  deck  represents  the  principal  engine 
simulation  model  used  by  the  Air  Force,  the  entire  deck  should  be  differ- 
entiated. 

It  is  difficult  to  estimate  cost  savings,  since  different  computer 
installations  use  different  cost  algorithms.  However,  it  is  possible  to 
estimate  processor  time  gains.  An  approximate  formula  for  estimating  pro- 
cessor time  gain  which  can  be  obtained  from  analytic  derivatives  was  derived 
during  this  program.  If  no  other  information  is  available,  processor  time 
gains  can  be  used  as  a rough  estimate  of  cost  gain.  The  derivation  of  this 
formula  together  with  several  examples  of  its  use  is  given  below. 

In  order  to  estimate  the  balanced  point  time  gain  due  to  the  inclusion 
of  analytical  derivatives  in  a gas  turbine  engine  cycle  deck,  the  following 
information  is  required: 

1.  The  number  of  base  points  per  cycle  balance  point. 

2.  The  number  of  derivative  points  per  cycle  balance  point. 

3.  Ratio  of  the  time  required  to  execute  an  analytical  derivavtive  pass 
to  the  time  required  for  a base  point. 

Items  1 and  2 can  be  obtained  from  the  iteration  history  of  the  deck 

averaged  over  a typical  flight  envelope.  The  third  item  may  be  estimated 

based  upon  past  experience.  For  example,  for  the  parametric  VCE  deck,  the 
ratio  is  known  to  be  about  one-tenth  (0.1). 

The  estimating  formula  is  obtained  as  follows: 

Let  NBP  “ number  of  base  points 

TBP  * processor  time  required  to  execute  one  base  point 

NDP  “ number  of  derivative  points  per  cycle  balance  point 

TOP  • processor  time  required  to  execute  one  analytic  derivative 
point 
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Assume  1.  The  time  required  to  execute  a base  point  is  10Z  greater  after 
analytic  derivative  logic  is  included  (due  largely  to  curve  and 
thermo  derivatives  which  must  be  evaluated  on  base  points). 

2.  The  time  required  to  execute  a base  point  is  the  same  as  the  time 
required  for  one  numerical  derivative  pass. 

Then, 

Time  for  Balanced  Pt . with  Analytic  Derivatives 
Time  Required  with  Numerical  Derivatives 


NBP  * TBP  * 1.1  + NDP  * TDP 
NBP  * TBP  + NDP  * TBP 

After  simplification, 


The  maximum  possible  time  savings  can  be  obtained  by  assuming  that  the  time 
for  an  analytical  derivative  pass  is  zero,  and  that  the  time  to  run  the  base 
point  is  unchanged  by  the  inclusion  of  the  analytical  derivative  code  (i.e., 
the  10Z  increase  in  base  point  run  time  does  not  exist).  With  these  two 
restrictions,  the  above  equation  reduced  to: 


Maximum  Time  Ratio  = 


1 


These  relationships  have  been  plotted  in  Figure  7.  The  percent  reduction  in 
run  time  in  the  present  deck  is  about  62Z  over  the  411  point  flighr  matrix. 
The  ratio  of  derivative  passes  to  base  points  averages  about  2.3.  The 
approximately  relationship  shown  in  Figure  7 would  predict  a reduction  of 
about  60Z  which  is  within  two  percentage  points  of  the  actual  value.  The 
maximum  time  gain  possible  would  be  70Z.  Therefore,  about  90Z  (62Z/70Z) 
of  the  possible  time  reduction  has  been  obtained  with  the  present  deck 
setup.  For  engines  and/or  flight  envelopes  having  a greater  ratio  of 
derivatives  passes  to  base  points,  considerably  greater  time  gains  can  be 
achieved  with  analytic  derivatives. 
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No.  of  Derivative  Passes/Vo.  of  Base  Points, 
(NDP/NBP) 


Figure  7.  Estimated  Reduction  in  Run  Time  Due  to  Analytic 
Derivatives. 
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For  iteration  algorithms  which  evaluate  derivatives  after  every  base 
point,  the  ratio  of  derivative  passes  to  base  point  passes  is  equal  to  the 
number  of  independent  variables.  For  the  VCE  computer  program  file,  the 
number  of  independent  variables  is  9.  Had  derivatives  been  taken  after  every 
point,  the  reduction  in  time  would  have  been  802  or  about  five  to  one. 

For  the  second  method  (low  level  derivative  set),  it  was  determined  that 
the  ratio  between  the  cime  for  a derivative  pass  to  a base  point  pass  is  0.29. 
The  curve  for  TDP/TBP  - 0.29  is  shown  on  Figure  7;  and  it  indicates  that  if 
derivatives  were  calculated  for  every  base  point  (NDP/NBP  ■ 9)  the  theoreti- 
cal gain  would  be  632.  These  results  compare  favorably  with  the  complete 
analytical  derivative  technique  and  the  method  requires  only  a minimal 
amount  of  time  and  effort  to  implement. 
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7.0  CONCLUSIONS 


It  can  be  concluded  that  in  any  engine  simulation  model  using  numerical 
derivatives  to  obtain  balanced  cycle  points,  a considerable  savings  in  both 
computer  time  and  cost  can  be  obtained  by  the  inclusion  of  analytic  deriva- 
tives. The  greatest  return  per  dollar  of  cost  results  when  only  that  portion 
of  the  model  for  which  the  coding  remains  relatively  constant  is  differenti- 
ated. For  the  Air  Force  engine  simulation  model  (SMOTE),  the  entire  deck 
should  be  differentiated  since  coding  changes  do  not  appear  to  occur  fre- 
quently . 

The  cost  savings  obtained  during  this  study  were  44%  of  a deck 
run  internally  and  52%  of  the  same  deck  ruu  externally  (the  differ- 
ence is  due  to  the  use  of  a larger  set  of  output  parameters  on  the  internal 
deck) . 
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